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Abstract 


The response of delaminated composite plates to compressive in-plane loads 
was investigated. The delaminated region may be either circular or elliptical, and 
may be located between any two plies of the laminate. For elliptical delaminations, 
the axes of the ellipse may be arbitrarily oriented with respect to the applied loads. 
A model was developed that describes the stresses, strains, and deformation of the 
sublaminate created by the delamination. The mathematical model is based on a 
two dimensional nonlinear plate theory that includes the effects of transverse shear 
deformation. The model takes into account thermal and moisture induced strains, 
transverse pressures acting on the sublaminate, and contact between the sublami- 
nate and plate. The solution technique used is the Ritz method. A computationally 
efficient computer implementation of the model was developed. The code can be 
used to predict the nonlinear load-strain behavior of the sublaminate including the 
buckling load, postbuckling behavior, and the onset of delamination growth. The 
accuracy of the code was evaluated by comparing the model results to benchmark 
analytical solutions. 

A series of experiments was conducted on Fiberite T300/976 graphite/epoxy 
laminates bonded to an aluminum honeycomb core forming a sandwich panel. Ei- 
ther circles or ellipses made from Teflon film were embedded in the laminates, sim- 
ulating the presence of a delamination. Each specimen was loaded in compression 
and the strain history of the sublaminate was recorded far into the postbuckling 
regime. The extent of delamination growth was evaluated by C-scan examination 
of each specimen. The experimental data were compared to code predictions. The 
code was found to describe the data with reasonable accuracy. 

A sensitivity study examined the relative importance of various material prop- 
erties, the delamination dimensions, the contact model, the transverse pressure 
differential, the critical strain energy release rate, and the relative growth direc- 
tion on the buckling load, the postbuckling behavior, and the growth load of the 
sublaminate. 
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Chapter 1 


Introduction 


Fiber-reinforced organic matrix composite materials may contain delamina- 
tions introduced, for example, by manufacturing defects or impact deimage. When 
subjected to compressive in-plane loads, the delaminated region may first buckle 
and then grow in size. Either of these occurrences may significantly limit the useful- 
ness of the composite plate. Therefore, to utilize the many inherent advantages of 
composite materials, the behavior of delaminations must be fully understood. This 
investigation addresses this problem, and specifically seeks to establish a model 
which predicts the buckling and postbuckling growth behavior of delaminations in 
composite plates subjected to in-plane compressive and shear loads. 

Owing to the significance of the problem, several investigators have proposed 
models describing the behavior of delaminated plates under compressive loading. 
The buckling and growth of through-width delaminations in plate strips have been 
analyzed by Chai et al. [1], Yin et al. (2), Simitses et al. [3], Gillespie and Pipes 
[4], Wang et al. [5,6], Sallam and Simitses [7], Williams et al. [8], El-Senussi and 
Webber [9], Vizzini and Lagace [10], Yin [11], and Kardomateas [12]. The buckling 
and growth of circular delaminations in isotropic plates under radial loads have 
been investigated by Bottega and Maewal [13], Yin and Fei [14], and Bruno [15]. 

In addition to these simple geometries, the behavior of plates containing rec- 
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tangular delaminations have been analyzed by Konishi [16] and Jones et al. [17], 
while plates containing elliptical delaminations have been investigated by Konishi 
[16], Chai and Babcock [18], Kassapoglou [19], Shivakumar and Whitcomb [20], and 
Whitcomb [21, 22]. In all but one of these analyses, the major axes of the rectangle 
or ellipse were assumed to be aligned with the direction of the compressive load. 
The one exception is the analysis of Shivakumar and Whitcomb [20], which assumed 
an arbitrary orientation for an elliptical delamination with respect to the applied 
load. However, Shivakumar and Whitcomb calculated only buckling loads and did 
not consider postbuckling deformation or growth. 

It appears that no model exists for predicting the buckling, postbuckling be- 
havior, and growth of: (a) circular delaminations or elliptical delaminations with 
axes arbitrarily oriented with respect to the applied in-plane loads, and (b) de- 
laminations located between any two plies of the composite plate. Therefore, the 
primary goal of this investigation was to develop a model capable of addressing this 
more general problem of delamination in a composite plate. 

In developing the model, an additional objective was to make the model readily 
useable in engineering practice. Substantial effort was made to keep the analysis 
simple and straightforward, to develop a computer implementation of the model 
that was computationally efficient, and to incorporate a user-friendly interface for 
the program so that meaningful results could readily be obtained. 

Finally, to firmly establish the credibility of the model, a series of experiments 
on graphite-epoxy face sheeted aluminum honeycomb sandwich panels containing 
embedded delaminations was conducted. Strain histories were measured far into the 
postbuckling regime, and the data were compared to the analytical results, verifying 
the validity of the model. 



Chapter 2 


Problem Statement 


Consider a multilayer laminated composite plate. The plate may be a “solid” 
laminate or a “sandwich” laminate consisting of two face sheets bonded to a honey- 
comb core (Figure 2-1). In either case, the plate must be symmetrically laminated 
with respect to its midplane. Each layer, or ply, in the plate may be made from 
a different material. Each material may be isotropic or orthotropic, the latter in- 
cluding continuous fiber reinforced composites. Each material must behave in a 
linearly elastic manner. A delamination exists between two adjacent plies in the 
plate (or face sheet) interior. The delamination may occur between any two plies, 
dividing the plate locally into two parts. The delamination may also be at the face 
sheet-honeycomb interface. Note that while the plate is symmetric, the two parts on 
either side of the delamination will, in general, be unsymmetric. The delamination 
is small with respect to the plate planar dimensions but large with respect to the 
thickness of the plate. The delamination is either circular or elliptical. The ellipse 
may have an arbitrary orientation with respect to the plate (Figure 2-2). 

In-plane tensile, compressive, and shear loads may act on the plate. The re- 
sponse of the plate to the applied loads is assumed to be initially unaffected by 
the behavior of the sublaminate formed by the delamination. The plate response 
to the applied loads determines the displacement boundary conditions for the sub- 


3 


Chapter 2: Problem Statement 


4 


laminate. Under the action of the applied loads, the delaminated sublaminate may 
buckle and subsequently grow in area. Given the plate material properties and the 
delamination geometry, the problem is to find: (a) the load-strain behavior of the 
sublaminate, (b) the load applied to the plate at which the sublaminate buckles, 
and (c) the load applied to the plate which causes an onset of delamination growth. 


1 1 
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Composite Plate Honeycomb Sandwich Plate 


Figure 2-1 Plate geometries investigated and the division of the plate into two 
parts as a result of the delamination. 




Chapter 3 


Delamination Analysis 


§3.1 Approach 

The major concepts of the analysis of a delaminated composite plate are pre- 
sented in this chapter. The analysis proceeds in four major steps. First, the displace- 
ments, strains, and stresses in the plate are calculated as though the delamination 
were not present. Second, the load at which the delaminated sublaminate buckles is 
determined. Third, the displacements, strains, and stresses in the sublaminate are 
determined using the condition that the displacements at the delamination bound- 
ary match those of the plate determined in step one. Fourth, the load at which the 
sublaminate grows in size is established. In the first step, the behavior of the plate 
is calculated directly from laminated plate theory [23]. In the remaining steps, the 
approximate behavior of the sublaminate is determined using energy methods. 

The following fundamental assumptions of plate theory are employed in the 
analysis: 

1. The thicknesses of both the plate and the sublaminate are small compared to 
all other dimensions. 

2. The thicknesses of both the plate and sublaminate are constant. 

3. The material behaves in a linearly elastic manner. 
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4. Each layer is either isotropic or orthotropic. 

5. The plate and the sublaminate undergo small strains, and the sublaminate 
experiences moderate rotations. 

6. The transverse normal stresses are zero in both the plate and the sublaminate. 

7. Perfect bonding exists between adjacent layers of the composite (except, of 
course, at the location of the delamination). 

8. The transverse displacements and rotations in the plate are zero. 

§3.2 Coordinate Systems 

Three cartesian coordinate systems are employed in the analysis (Figure 3- 
1). The coordinate system coincident with the principal material axes of each 
orthotropic ply is the x, y, z system. The coordinate system coincident with the 
semi-axes of the delamination ellipse is the x \ , r 2 , x 3 system. The coordinate 
system of the plate is the xy, xy, xy system. The r, y, z system is the on-axis 
system, while the other two systems are off-axis systems. The two off-axis systems 
are related to the on-axis system by rotations about the transverse axes, where 
these axes are all equivalent (z = x 3 = xy ). 

§3.3 Constitutive Relations 

The constitutive relations for a linearly elastic material are 

(Tij = Cijkitki ij = 1,2,3 (3.1) 

where the Cijui are elastic constants relating the stresses <7,j to the strains e*/. For an 
orthotropic or isotropic material in the on- axis coordinate system, the constitutive 
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Sublaminate 


Figure 3-1 Cartesian coordinate systems: xyxyxy for the plate, X1X2X3 for the 
sublaminate, and xyz for the material coordinates of each ply. 6 is 
the angle between the plate and sublaminate systems, and <f> is the 
angle between the plate and ply systems. 
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(3.2) 

where the ot{ are the on-axis linear coefficients of thermal expansion. The temper- 
ature difference AT is defined as 


A T = T-T r 


(3.3) 


where T is the uniform temperature of the composite and T r is a reference temper- 
ature at which the thermal strain is defined to be zero. A convenient value for T r 
is the temperature at which the material “solidifies” during curing. 

The effects of moisture absorption by the composite material can be treated 
in an analogous manner. The strain due to moisture uptake is /?, Ac, where the fli 
are the on-sods linear coefficients of moisture-induced expansion and Ac is the rel- 
ative change in moisture concentration. For simplicity, the analytical development 
presented here is in terms of thermal strsuns. However, an equivalent analysis of 
moisture effects can be made by substituting /?,Ac for OjAT. 

Contracted notation will be used in the rest of the analysis except where noted. 
For example, in an off-sods x x , x 2 , and x 3 coordinate system the stresses, strains, 
sind elsistic constsmts are represented by [24] 
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Note that e 4 , e 5 , and are engineering strains. 

For a material in plane stress, the constitutive relations (Eq. 3.2) may be 
conveniently written in an off- axis coordinate system as 
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where the Qij are the plane stress reduced stiffnesses defined as 

Qij = C„ - ij = 1,2,6 

^33 

Qli = C„ i,j= 4,5 


(3.5) 


(3.6) 


The stresses and strains have been arranged so as to group the in-plane and out- 
of-plane components separately. The apparent thermal shear term (at AT) appears 
since the constitutive relations are expressed in an off-axis coordinate system. Ex- 
pressions for Qij in terms of engineering constants for isotropic and orthotropic 
materials are given in Appendix A. 


§3.4 Displacements, Strains, and Stresses in the Plate 

In the plate coordinate system (x\i,xy,xy), the in-plane total strains in the 
symmetrically laminated plate (containing no delamination) resulting from uni- 
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fonnly applied in-plane load resultants Ny, N 2 >, and N 6 > (Figure 3-2), and a tem- 
perature difference AT are [23] 

/ \ / dyy Oy 2 > a l'6' \ / Ny + Nf, \ 

°4 I = a V2 . a vv a 2 . 6 , N 2 . + iV|/ T (3.7) 

\°e%!/ \ a l'«' a W ae'6'/ \Nv+N’! t J 

The matrix elements ay, which axe the inverses of the plate stiffness matrix elements 
Aij [23], axe 


(«»'>') = 1 = 



Qi'j' 



i'J' = 1,2,6 


(3.8) 


where h pl is the thickness of the plate. The thermal load resultants Nf! T , N*} T , 
and Nl ,T are 



Qi> 2 > 
Q vv 
Qw 


Qvv 
Q 2 ' 6 ' 
Q&' 6' 



AT<fx3# 


(3.9) 


In the sublaminate coordinate system (xi, x 2 , X3), the total strain components 
are determined by a tensor rotation about the transverse axis X 3 > 

/°<\ / ™ 2 n 2 nm \ /°eff\ 

( •ej 1 1 = ( n 2 m 2 nm I I °c%! I (3.10) 

\«eJV \-2nm 2nm m 2 - n 2 / \°cjf/ 

where m and n are cos 6 and sin#, respectively, and Q is the angle between the plate 
and sublaminate coordinate systems (Figure 3-1). 

To simplify the buckling and postbuckling analyses, proportional mechanical 
loading is assumed. Each of the in-plane load resultants may have a unique value, 
but the relationship of one to another is fixed. In this way, a single load parameter 
N suffices to characterize the total load on the plate 
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N v \ / 7l 

N v = 72 | N 

■^ 6 '/ \ 76 , 


( 3 . 11 ) 


with the relative magnitude of each load described by the fractions 71 , 72 , and 76 
(71 +72 + 76 = 1)- The resulting total strains in the sublaminate coordinate system 
(xj, X 21 £3) may now be expressed as 



N + 


' o„plT > 
€ 1 

E PI 

h_ 

o,P 

e 6 


\ PIT 


IT 


where cj, C2, and C3 are defined by 



m‘ 


n* 


nm 

n“ m* nm 
—2 nm 2 nm m 2 — n 2 


and the thermal strains in the plate by 



m* 


n* 


nm 

2 2 

n m* nm 

—2 nm 2nm m 2 — n 2 


aw 

ai»2' 

<*1'6' 

a 1*2' 

a 2 * 2 * 

aw 

<*1'6' 

ay# 

«6'6' 


ayy 

a 1*6' 

ai>2> 

<12*2* 

a w 


aw 

a 6 / 8 # 





(3.12) 


(3.13) 


Integration of Equation 3.12 (Appendix B) gives the in-plane displacements of the 
plate 


X' = (ciJV + •«5* r )*i + i(caJV + 
x' = l(ceJ V + “<; ,T )x, + (c 2 W + 


(3.15) 


Equation 3.15 together with Assumption 8 (°u£ 7 = 0) completely describe the dis- 
placements in the plate. The displacements at the sublaminate boundary must 
match these displacements. Note again that the displacements and strains are 
taken to be zero when the temperature is T r and no mechanical loads act on the 
plate. 
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The nonlinear strain-displacement relations used for the sublaminate are those 
proposed by von Karman [25] for the large displacement analysis of plates. Although 
von Karman discussed the use of these strain-displacement relations only in the 
context of classical plate theory, it has been shown (Appendix C) that they are 
appropriate for the moderate rotation, shear deformation theory used here. The 
strains (reverting to conventional notation for the moment) are related to the 
displacements by 


— Qt* 4. £lfi i d “3 du 3 , 

2 dij dn dx{ dxj 


(3.16) 


Using a higher order shear deformation theory [26-29], the sublaminate dis- 
placements are taken to be cubic functions of the transverse coordinate x 3 


(* 1 , * 2 , x 3 ) = °uf (xi, x 2 ) + x 3 ipi{xi,x 2 ) + , x 2 ) + x\r)[\x x , x 2 ) 

7 x 2 7*3) = °«2 (Xi,X 2 ) -I- X^\x u x 2 ) + x|^ , (x 1 ,x 2 ) -I- x^^(x x ,x 2 ) 

*f(*u*2 t *») = °^z{x\,x 2 ) 

(3.17) 

where “uj*, 0 u^, and °u\ l are the displacements of the sublaminate midplane, and 
V’f , 7 £ 1*7 'lit *12 describe rotations of line segments originally perpendicular 

to the sublaminate midplane. 

No shear forces act on any of the lateral surfaces of the sublaminate. Hence, 
the shear stress components (returning to contracted notation) on these surfaces 
are 

h * 1 

<74 = as = 0 at x 3 = ±— (3.18) 


where h al is the thickness of the sublaminate. For a sublaminate constructed of iso- 
tropic or orthotropic materials, the above condition requires that the corresponding 
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shear strains on the lateral surfaces also be zero (Eq. 3.2). 


€4 = 65 = 0 at 



(3.19) 


Substituting the displacements (Eq. 3.17) into the strain-displacement relations 
(Eq. 3.16), differentiating, and applying the four boundary conditions (Eq. 3.19), 
the relationships between the rotation functions are [30] 


(!' = «' = 0 



3 (k“) 

4 


, dui 1 , 


Z(h‘‘) 




(3.20) 


Using Eq. 3.20, the displacements (Eq. 3.17) may be rewritten as 


= •«;' + *3 W 


4xl , du\ l 


3 (h* l Y y dx 






• ill — 
u 3 = U 3 


(3.21) 


The displacements are now specified in terms of only five functions: three midplane 
displacements °u* 1 , *^3*, and two rotation functions tp * 1 and ij >2 • Using the 

displacements in Eq. 3.21, the nonlinear strains (Eq. 3.16) may be expressed as 


e[ l = °e[ l + *3 + x 3 3 2 k { 1 

e 2 = ° e 2 + * 3 0 « 2 * + xl 2 K2 
t\ l = 0 

(3.22) 

4 = *4 + A 2 «;' 

4 = °4 + A 2 4 
4 = *4 + x 3°4 + *3 *4 
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In these expressions, °e * 1 are the midplane strains, and and are the midplane 
curvatures of the sublaminate defined by 


,f fjV . 1 ( 

d Xl 


«i = 


d°u\ l 


dx\ 




vi ax' i,9v' 
2 &T2 2^ dx 2 


3 )2 




= v-r + 


o-*/ _ 


= i>l l + 


d%l 

dx 2 

dxi 


V ,_9X' 


9X' 


+ 


5xi 
d°u^ d°u\ l 
dx\ di 2 


.1 w 

dxi 

diii' 


*i - 


•*5'- 


X 

5^2 


« 6 — 


dl 2 


+ 


dx\ 


2 _ 
K \ — 


2 al _ 
rx 2 — 


2„W _ 
« 4 — 


2*1 _ 
« 5 — 


2 „»l _ 

«« — 


3(^0 

4 


.0 2 °^' , drfil 


dx x ‘ 


+ — —) 
+ dx x ’ 


3(^0 

4 


,d 2 o ui l 

~2 ( «■» 9 + 


di 2 2 


ur 50^ 


5x^ 


* . ^ i 3 < \ 

(S?'* + W 
(I=F ( * + ‘KT* 

4 g< djff 
) 


2(h tl y y dx 2 


+ 2 


a 2o u| f 

dx\di 2 


(3.23) 


§3.6 Displacements, Strains, and Stresses in the Sublaminate 

The displacements in the sublaminate are assumed to be a linear combination 
of two parts: (a) the displacements that would exist in the sublaminate in the 
absence of the delamination, plus (b) the displacements introduced by transverse 
deformation of the sublaminate. The five functions describing the sublaminate 
displacements are 


X' — X' + ”4;' 

X' a x' + x' 

X' a X' + X' = X' (3.24) 

< ax' + «' = «' 
x' = x' + X' = X' 
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The in-plane displacements in the plate (°Uj , 1 °u% 1 ) are given in Eq. 3.15. The second 
equalities in the last three expressions can be written because the transverse dis- 
placements and rotations in the plate Eire assumed to be zero (Assumption 8). The 
displacements in the sublaminate due to out-of- plane deformation (the quantities 
with the hat) are represented by 


j = 1 

>=1 

"“ 5 1 

y= i 

ft* = £•/*< 
*2 = 

j = i 


(3.25) 


where n p through n t are the number of terms in each series. The parameters pj, qj, 
rj, Sj , and tj are coefficients, while p <f>j, q 4>j, r 4>j, and *<f>j are fimctions of the 
coordinates x\ and x<i- Expressions for these coordinate fimctions must be chosen 
such that they: (a) satisfy the boundary conditions on the sublaminate (discussed 
below) and (b) are linearly independent, continuous, and complete [30]. 

One of the f undam ental assumptions of this analysis is that the displacements of 
the sublaminate and plate match along the boundary of the sublaminate. Along this 
boundary, the displacements are completely specified, while the force and moment 
resultants are unspecified. Thus, by definition, the sublaminate boundary is a 
clamped bo un dary (Figure 3-3). The appropriate clamped boundary conditions are 
[26, 27] 
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Clamped Boundary 



Figure 3-3 Illustration of the clamped boundary at the sublaminate edge. 


• u ;'="u;' = o 


«' = «' = o 


v>;' = v-f' = o 

a°n;' 

9x n 


a°u?' A 

dx n 


(3.26) 


The subscripts n and t refer to the normal and tangential directions, respectively, to 
the sublaminate boundary (Figure 3-4). Clearly, the first five boundary conditions 
require that the functions p <f>j, q <f>j, r <f>j, and *4>j vanish at every point on the 
boundary (Eqs. 3.24 and 3.25). The final boundary condition of Eq. 3.26 requires 
that the derivative of r <f>j also vanish on the boundary. Accordingly, the following 
polynomial coordinate functions were chosen 


Chapter 3: Delamination Analysis 


19 



Figure 3-4 Definition of the tangent t and normal n coordinates along the sub- 
laminate boundary. 
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= (1“*1 -*a)(Pl*l + P 2*2 + P3*? + P4*?*2 + P5*1*2 + P 6*2 
+ P7*i + P 8*?*2 + P9*l*2 + PlO*i *2 + Pll*l*2 + Pl 2*2 
+ Pl3*J + Pl4*i*2 + Pl5*i*2 + Pl 6*?*2 
+ p 17 XiX^ +PI 8 X 1 X* + Pl9*l*2 +P20^2) 

qj 9 <t>j = (1 - x\ - + 92*2 + 93*1 + 94 * 1*2 + 95 * 1*2 + 96*2 

+ 97*1 + 98*1*2 + 99*1*2 + 910*1*2 + 911*1*2 + 912*2 

(3.27) 

+ 913*1 + 914*1*2 + 915*1*2 + 916*1*2 
+ 917*1*2 + 918*1*2 + 919*1*2 + 920 *2) 


r j r< t>j = (1 “ *1 “ *!) 2 ( r l + 1 * 2*1 + 1*3*2 + r4*lZ2) 


Sj* 4 >j — (1 — x\ — x|)(3iXi + 32*2 + 83X* + 34X1X2 + 35X1X3 + 36x|) 


= (1 — Xi — x|)(*ixi + < 2*2 + < 3*1 + < 4 X 1 X 2 + < 5 X 1 X 2 + < 8 * 2 ) 


where 


r x i\ 
/ *2 

x, = ( T ) 


These expressions satisfy the boundary conditions for an elliptically-shaped sublam- 
inate (semi-axes a and 6), and satisfy the condition that the functions be linearly 
independent, continuous, and complete. 

The functions in Eq. 3.27 are similar to those used by previous investigators 
[16,18,19]. However, previous investigators have omitted various terms from the 
series. Either the omission of terms from a series or the premature truncation of a 
series can affect the accuracy of the results. In particular, it is important to retain 
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1. Crossproduct terms (e.g. 74 X 1 X 2 ) for arbitrary delamination orientations; 

2. p 4>j and q <f>j to at least one order higher than r <f>j for accurate postbuckling 
calculation of in-plane strains; and 

3. *<f>j and *<f>j to the same order as for accurate representation of transverse 
shear rotations. 

The displacements in the sublaminate (Eq. 3.21) are now specified in terms 
of the unknown coefficients pj, qj , ry, Sj, and tj and their associated coordinate 
functions (Eq. 3.27). The off-axis strains at any point in the sublaminate may 
be calculated from the midplane strain and curvature definitions (Eq. 3.23) and 
the strain-displacement relations (Eq. 3.22). Using the strains, the associated 
off-axis stresses in each ply of the sublaminate may be calculated from the off-axis 
constitutive relations (Eq. 3.5). On-axis strains and stresses in each ply are obtained 
by rotating the respective off-axis strains and stresses into the on-axis coordinate 
system [23]. 

The displacements, strains, and stresses in the sublaminate are thus specified 
in terms of the coefficients pj, qj , rj, sj , and tj. These coefficients are determined 
by the Ritz energy method. 


§3.7 Total Potential Energy 


The total potential energy of the sublaminate n* < in the absence of body forces 
is [31] 




n 4 ' = jQ (Tide^jdV - j(^J fdu^JdA i = 1,2, 4, 5, 6 (3.28) 


VO AO 

where V is the volume of the sublaminate, A is the lateral surface area of the 
sublaminate, and / is the force per unit area acting on the surface. The subscript 
3 is not included since the transverse normal stress <73 is assumed to be zero. Note 
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that repeated subscripts imply summation. Substituting the constitutive relations 
(Eq. 3.5) for <J{ yields 



„« i 


ajAT)de^ dV - jQ f du^ dA 


ij = 1,2, 4, 5, 6 (3.29) 


Integration of Eq. 3.29 with respect to e; produces 


n -1 = 


\ J eiQaejdV- J CiQijatjAT dV 

V V 


J(ff du3 ') dA 

A 0 


1,2, 4, 5, 6 

(3.30) 


Two kinds of surface traction are considered to act on the lateral surfaces of 
the sublaminate. First, a uniform transverse pressure may exist due to a pressure 
difference A P between the outside and inside surfaces of the sublaminate (Figure 3- 
5). The pressure on the outside surface is generally atmospheric. On the inside 
surface the pressure may be subatmospheric due to a partial vacuum that may 
form as the sublaminate buckles. Second, a force may result from contact between 
the buckled sublaminate and the plate over portions of the delaminated surface 
(Figure 3-6) [32]. Where the sublaminate tends to deform toward the plate, contact 
between the two will occur and a force resisting the sublaminate deformation will 
arise. This contact is modeled by considering the sublaminate to be resting on a 
detached elastic foundation (Figure 3-7). The restoring force is taken to vary linearly 
with the sublaminate transverse displacement uj* for positive displacements, and to 
vanish for negative displacements. Therefore, the force per unit area acting on the 
sublaminate at a given point is 


, / AP - Kul l ui l > 0 

/ = {ap 4'<o 


(3.31) 
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The foundation modulus K [33, 34] is estimated from (Appendix D) 


K 



(3.32) 


where Ef is the elastic modulus of the foundation and If is a characteristic length 
(for example, the sublaminate diameter). 

Note that the transverse pressure model only makes sense in conjunction with 
the contact model. Without the contact model, the transverse pressure would 
simply produce sublaminate bending toward the plate. 


Atmospheric Pressure 



Figure 3-5 Possible pressure difference acting across the sublaminate thickness. 
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Figure 3-6 Contact between the sublaminate and plate. 
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Integration of Eq. 3.30 with respect to the thickness h tl , together with the 
strain definitions (Eqs. 3.22 and 3.23), results in the following expression for the 
total potential energy [30] 


n" = 



A 45 £44 £45 

A 55 D\s D 55 
sym F 44 F 4S 
F 55 
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/ n; ,t \ 

/ n *it \ 

Nj lT 

M{ lT 

Mj lT 
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p) lT 

V Pi ,T / 


dA 


*•1 


-JU /du ’) 

A 0 


dA 


(3.33) 


where tr represents the matrix transpose, and A,y, Z?,y, Z?,y, 2?,y, 2^-y, fTy are 
sublaminate stiffnesses (Appendix E) defined by 
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* 


(,Aij i Dij ) Eij f Fij , J X3, X3, £3, £3, £3) dxj i^j — 1 , 2,6 




*r 

(Aij, Dij, Fij) = J Qij(\,x\,x\)dxz i,j = 4,5 


-Hr 


The thermal force and moment resultants are 




(3.34) 


(IV * ,r , M* lT , P* lT ) = J Q ijaj AT(l,x 3 ,xl)dx 3 ij- 1,2,0 (3.35) 

The specific limits of integration over the area A of the elliptical sublaminate 


are 


r* fby/l-W 
/ / __[• • •] dx 2 dxi 

J —a J —b\/ 1 — (-?) 


(3.36) 


r« [by/ 

* 2 2 

where the ellipse is bounded by 1 — ( ^-) — (^) =0, and a and b are the semi-axes 
in the x\ and x 2 directions, respectively. 


§3.8 Applied Load versus Deformation of the Sublaminate 

To establish the relationship between the applied mechanical and thermal loads 
(N and AT) and the sublaminate deformation (ti**, u^, and uj 1 ), the total potential 
energy of the sublaminate (II'*) is first assembled by substituting the displacement 
approximations (Eqs. 3.24, 3.25, and 3.27) into the midplane strain and curvature 
definitions (Eq. 3.23), and then substituting these results into the expressions for the 
total potential energy (Eq. 3.28). The resulting expression is extremely lengthy and 
will not be given here. Essentially, the total potential energy of the sublaminate 
is now expressed as a function of the known geometry (a, 6), material properties 
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(Aij , . . . , a,), and applied loads (N, AT), and a set of as yet unknown coefficients 

(Pji 9;, rj, Sj, tj) 


II — II (a,,b] Ajj, . . . ,an‘, N, &T',pj,qj,r j, 


(3.37) 


The unknown coefficients axe determined by minimizing the total potential energy 
with respect to the coefficients [30] 


dll** 

drrij 


= 0 j — 1 to 56 


(3.38) 


where my is used as a generic unknown coefficient representing py, ?y, ry, sy, or tj. 
The differentiations indicated in Eq. 3.38 result in a system of 56 nonlinear alge- 
braic equations in the unknown coefficients. These are the equilibrium equations. A 
solution to these equations yields a set of coefficients (py,qy,ry,iy,ty) corresponding 

A A 

to specified values of the applied loads (N, AT), where the hat indicates a partic- 
ular set of loads. Knowing the coefficients, the displacements, strains and stresses 
throughout the sublaminate can be determined. 

A load-deformation history for the sublaminate is mapped out by solving the 
equilibrium equations over a range of applied loads. However, care must be exercised 
due to the nonlinearity of the equations. In general, more than one solution exists for 
a given load. The solution must correspond to a local minimum of the total potential 
energy, implying that the solution must be physically stable. Furthermore, multiple 
stable solutions are possible. Therefore, each possible solution (py, 9y, fy, 5y,ty) must 
be tested to determine whether it corresponds to a local minimum of the the total 
potential energy (stable solution) or a local maximum (unstable solution). The 
stability test [35] requires that the determinant of the matrix of second partial 
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derivatives of the the total potential energy be positive definite 


d 2 n**(pj, 


dmidrrik 


> 0 i, k = lto56 


(3.39) 


where m,- is again a generic coefficient. If multiple stable solutions axe shown to 
exist, each with an associated total potential energy, then the solution with the 
minimum total potential energy of the sublaminate is assumed to be the most likely 
solution. 


§3.9 Buckling Condition 


The buckling load is one at which the system changes from one configuration 
to another, energetically more favorable, configuration. This change occurs at a 
load for which the determinant of the matrix of second partial derivatives of the 
potential energy ceases to be positive definite [31, 36] 


drriidmk 


= 0 


i, k — 1 to 56 


(3.40) 


In the equilibrium problem described above, the unknown coefficients are deter- 
mined for a known load. In the buckling problem, both the load and the coefficients 
for which Eq. 3.40 applies are desired. Two different approaches have been used 
to solve the buckling problem. In the linear method, the values of the unknown 
coefficients are assumed to be zero (pj = q 3 = rj = Sj = tj = 0) and the buckling 
load Ny satisfying Eq. 3.40 is found. In the equilibrium method, the load N is 
gradually incremented over a range of values, each time solving the nonlinear equi- 
librium equations for the unknown coefficents as above. The point at which the 
displacements change dramatically with increasing load (Figure 3-8) is the buckling 
load Nh. Since unsymmetric laminates may deform out-of-plane at loads less than 
the linear buckling load, the equilibrium method is recommended. 
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Figure 3-8 Illustration of the actual load-strain behavior and the calculated linear 
buckling load. 


§3.10 Growth Criterion 

The strain energy released per unit area by the plate-sublaminate system for 
an increment of sublaminate growth is the strain energy release rate G. The delam- 
inated sublaminate is assumed to grow for a given load when G exceeds a critical 
strain energy release rate G c of the material [37] 


Cs<i (ir'-n") 

dA 


(3.41) 


where IP* and IP* refer to the strain energies of the sublaminate and plate, respec- 
tively, and A is the surface area of the sublaminate. The strain energy released by 
the plate is the strain energy of that portion of the plate which becomes part of the 
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new sublaminate after growth (Figure 3-9). Thus, the strain energy of the plate is 


ip' 


-\Si I 



A\2 A.22 A.26 

•^16 ^26 ^66 


o,pl 


dA 


- 1 dA 


(3.42) 


N?r 


The stiffnesses A{j and thermal resultants Nf ,T are evaluated over the thickness 
of the sublaminate because the strain energy of the balance of the plate does not 
change during growth of the sublaminate. Only the sublaminate portion of the 
plate contributes. 

Following Chai and Babcock [18], the total strain energy release is considered 
during growth of the sublaminate. That is, although the strain energy released 
during growth of the sublaminate varies along the sublaminate boundary, local 
variations in the strain energy release are not included here. For an elliptically- 
shaped sublaminate, Eq. 3.41 may be written (Appendix F) 

■ an* 1 da ■ an* 1 \ /air* da , air* 


/dll'’ da 1 dll'’ \ /dll” da 1 air' \ 

Q = (sriE + ST) - (-5T7E + ST) > Q 

+ <■) 


(3.43) 


where A = irab is the area of the elliptical delamination. The parameter describes 
the direction in which the sublaminate grows (Figure 3-9). For example, = 0 
implies growth in the b direction only, ^ — 00 implies growth in the a direction only, 
and j implies self-similar growth. The strain energy release rate is evaluated 
over a range of values of the parameter ^ so that the lowest G can be found. In 
practice, sublaminate growth is often observed in a direction perpendicular to the 
applied load. Thus, a suitable choice would be ^ = 0 (growth in the b, or X 2 , 
direction) for a load applied in the x x direction. 
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Chapter 4 


Implementation 


§4.1 Introduction 

The FORTRAN computer code DELAM was developed from the analysis of 
sublaminate buckling and postbuckling behavior (Chapter 3). The program reads 
input data describing the delamination and the plate in which it is contained, in- 
cluding: (a) the plate material properties, geometry, and layup; (b) the location, 
dimensions and orientation of the delamination; and (c) the applied loads. A list of 
the required input data is given in Table 4-1. 

From the input data, derived properties are calculated for subsequent use in 
the delamination analysis: (a) the on-axis stiffnesses of each ply (Appendix A); 
(b) the ply on-axis plane stress reduced stiffnesses (Eq. 3.6); (c) the ply off-axis 
reduced stiffnesses which appear in Eq. 3.5 in both the plate and the sublaminate 
coordinate systems; (d) the laminate stiffnesses for the plate (Eq. 3.8) and the sub- 
laminate (Eq. 3.34) from the off-axis reduced stiffnesses; (e) the thermal resultants 
for the plate (Eq. 3.9) and sublaminate (Eq. 3.35) from the off-axis reduced stiff- 
nesses, the thermal coefficients of expansion (as rotated into the off-axis system 
[24]), and the specified temperature difference; (f) the constants describing the me- 
chanical response of the plate (Eq. 3.13); and (g) the thermal strains (Eq. 3.14). 
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Table 4-1 Input Parameters Required for the DEL AM Computer Code 


Plate and sublaminate geometry and layup 
Number of plies in the plate, k p 
Number of plies in the sublaminate, 

Semi-axes of the ellipse, a, b 

Angle of the sublaminate axes with respect to the plate, 9 
Thickness of each ply, ti 
Orientation of each ply, 4>\ 

Material properties for each ply 

Longitudinal Young’s modulus, E x 
Transverse Young’s modulus, E y 
Longitudinal to transverse Poisson’s ratio, u zy 
In-plane shear modulus, G xy 
Out-of-plane shear moduli, G tz ,G yx 

Thermal (or hygro) coefficients of expansion, a x ,a y (f3 x ,/3 y ) 

Growth and contact parameters 

Critical strain energy release per unit area, G c 
Relative growth direction parameter, 

Contact law foundation modulus, K 
Load description 

Normal load in the 1 direction, 71 
Normal load in the 2 direction, 72 
Shear load in the 1-2 plane, 76 

Change from reference temperature (or from dry) state, AT (or Ac) 
Transverse pressure load, A P 
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The total potential energy of the sublaminate is calculated from Eq. 3.33, using 
the laminate stiffnesses and thermal resultants which have already been determined. 
The substitutions and integration were performed by the symbolic mathematics 
program MACSYMA [38], with the exception of the contact model, which was in- 
tegrated numerically [39]. An expression for the total potential energy was thus 
established as a function of the known geometry, material properties, applied loads, 
and a set of unknown coefficients to be determined. The first and second mixed 
partial derivatives of the total potential energy with respect to the unknown coeffi- 
cients were determined using MACSYMA. In addition, the partial derivatives of the 
total potential energy with respect to the major and minor axes of the sublaminate 
ellipse were evaluated using MACSYMA. The expressions for these derivatives were 
then inserted into the computer code DELAM for use in the load-strain behavior, 
buckling, and growth calculations. 

§4.3 Nonlinear Load-Strain Behavior 

The stresses and strains in the sublaminate are determined by obtaining solu- 
tions to the equilibrium equations (Eq. 3.38) for specified values of the load N. The 
equilibrium equations are a set of 56 simultaneous algebraic equations nonlinear 
in the coefficients m,-. These equations are solved for the unknown coefficents m, 
by the Newton- Raphson method [39]. Once a solution for a given N is found, the 
stability of the solution is checked by calculating the determinant of the matrix of 
second partial derivatives of the total potential strain energy using the values of the 
coefficients obtained in the solution (Eq. 3.39). A positive determinant indicates 
that the solution is stable. 

The displacements at any point in the sublaminate are calculated by substi- 
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tuting into Eq. 3.24 the solution coefficients and associated coordinate functions 
(Eqs. 3.25 and 3.27) together with the boundary conditions (Eq. 3.15). The mid- 
plane strains and curvatures are then determined from Eq. 3.23. The strains at any 
point in the sublaminate are calculated from Eq. 3.22. Finally, the stresses associ- 
ated with these strains are determined from the constitutive relationship (Eq. 3.5). 
The load versus strain behavior is determined by repeating the above procedure for 
different values of the applied load N. 

§4.4 Buckling Load 

The linear buckling load N l b is the load at which the determinant of the matrix 
of second partial derivatives of the total potential energy equals zero (Eq. 3.40). 
Using the given geometry and material properties, the elements of the matrix are 
numerically evaluated for an initial estimate of the buckling load, assuming that the 
unknown coefficients are equal to zero. The determinant of the matrix is found by 
decomposing the matrix into lower and upper triangular matrices. The product of 
the diagonal elements of the upper triangular matrix (LU decomposition [39]) is the 
value of the determinant . This constitutes a single evaluation of the determinant as 
a function of N. In general, the determinant is a nonlinear function of the load N, 
and explicit derivatives of the function with respect to N do not exist. The load at 
which the determinant is zero is found using the secant method [39]. 

Alternatively, the buckling load Nb is graphically determined by examining the 
complete load-strain behavior of the sublaminate. The load at which the sublam- 
inate behavior begins to markedly deviate from a linear response is defined to be 
the buckling load. 



Chapter 4: Implementation 

§4.5 Growth Load 


37 


The growth load of the sublaminate is the load at which the strain energy 
release rate of the plate-sublaminate system exceeds the critical strain energy release 
rate of the material (Eq. 3.41). A value is assumed for the growth load N g , and the 
associated nonlinear equilibrium displacements are determined in the same manner 
as in Section 4.3 above. For the displacements thus obtained, the derivatives of the 
total potential energy with respect to the geometry are evaluated and the strain 
energy release rate G of the system is calculated (Eq. 3.43). This G is a nonlinear 
function of N. The value of N at which G equals G c is found using the secant 
method [39]. 

§4.6 Code 

The computer program DELAM was specifically written to be used for design 
calculations as well as for research. It has a user friendly interface, and is computa- 
tionally efficient and fast. For example, the computation of the nonlinear load-strain 
behavior of a sublaminate over sixteen values of the applied load requires 7 minutes 
of CPU time on a Sun 3/160 workstation. The input parameters required by the 
code are given in Table 4-1. The outputs provided by the code are listed in Table 4- 
2 and illustrated in detail in Chapter 8. Sample input and output of the code are 
included in Appendix G. 
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Analytical Verification 


§5.1 Introduction 

Verification of the delamination behavior model consists of three tasks: (a) com- 
parison to known benchmark analytical solutions, (b) comparison to other approxi- 
mate solutions, and (c) comparison to experimental data. The first two verification 
tasks are presented here; the experimental procedure and results are presented in 
Chapters 6 and 7, respectively. 

The benchmark problems consider the behavior of circular and elliptical plates 
without delaminations under various loadings. The computer program DELAM 
must be able to predict the behavior of simple plates under edge compression and 
uniform pressure loads. For certain geometries and material properties, closed form 
analytical solutions exist. These have been chosen as the benchmark problems. 

Approximate solutions to the behavior of plates containing elliptical delami- 
nations have been proposed by several investigators [16,18,19,20,21,22]. In general, 
these solutions pertain to problems more limited than the analysis presented here. 
Nevertheless, some of the approximate solutions may be compared to the present 
method for a select set of problems. 
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§5.2 Buckling of Circular and Elliptical Plates Without Delaminations 

Consider the buckling of an isotropic, circular plate subjected to uniform edge 
compression. The buckling coefficients k, defined as 

k = N b ^ (5.1) 

were calculated by DELAM for both clamped and simply supported aluminum 
plates, where N b is the critical buckling load and a is the plate radius. The plate 
bending stiffness D is 

Eh 3 

~ 12(1 - i/ 2 ) 

where E is Young’s modulus, h is the plate thickness, and v is Poisson’s ratio. The 
classical buckling coefficients for clamped and simply supported plates are given by 
Timoshenko and Gere [33] as 14.68 and 4.20, respectively. The coefficients are inde- 
pendent of the plate radius-to-thickness (aspect) ratio f since the solution is based 
on classical plate theory. The DELAM and classical buckling coefficients are plot- 
ted versus the plate thickness ratio in Figure 5-1. The primary difference between 
classical plate theory and the present method is the inclusion of transverse shear 
deformation in DELAM. At large thickness ratios, the predictions are identical; at 
thickness ratios of less than 20, the effects of shear deformation become apparent 
as the DELAM buckling coefficient drops significantly below the classical value. 

Consider next the buckling of an isotropic elliptical plate under uniform edge 
compression. The buckling coefficients k for clamped and simply supported plates 
were calculated by DELAM for an aluminum plate over a range of ellipticities (j) 
from one to five. The present results are compared with tin approximate solution 
developed by Voinovsky-Krieger [40] in Figure 5-2. As Voinovsky-Krieger did not 
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Figure 5-1 Buckling coefficients of clamped and simply supported circular alu- 
minum plates subjected to uniform edge compression (a = 1.0 in.). 



Figure 5-2 Buckling coefficients of clamped and simply supported elliptical alu- 
minum plates subjected to uniform edge compression (a = 1.0 in., 
h = 0.01 in.). 
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include shear deformation effects, a large thickness ratio ( f = 100) was chosen for 
the present calculation to minimize the effects of shear deformation. At the lower 
ellipticities the results are virtually identical, while at higher ellipticities the present 
results are slightly lower. This is to be expected since the present solution uses more 
terms in the approximating functions than did Voinovsky-Krieger, thereby reduc- 
ing the stiffness of the approximation and lowering the buckling coefficients. The 
simple support prediction is included in the figure for completeness since Voinovsky- 
Krieger suggested that the buckling coefficients for simply supported plates could 
be estimated by dividing the clamped plate values by a factor of 3.5. 


§5.3 Large Deflections of Circular Plates Without Delaminations 


The linear theory of plate bending is usually limited to transverse deflections 
on the order of fractions of the plate thickness. Nonlinear theories which include 
moderate rotations (such as the present method) allow transverse deflections up to 
about two times the plate thickness. 

Consider an isotropic, clamped circular plate without delaminations subjected 
to a uniform transverse load. The transverse deflections at the center of the plate 
were determined by DELAM as a function of the applied uniform load. In Figure 5- 
3, the present solution is compared to a perturbation method solution given by Chia 
[41] 


goQ 4 

Eh* 


16 


tv 0 


3(1 - 1/2) h 


1 


,w a 


1 + 35o (1 + ,/,(173 - 73 " )( T )2 


(5.3) 


where q 0 is the transverse uniform pressure and w 0 the transverse deflection at the 
center of a thin plate. The present solution and Chia’s solution agree very well 
out to transverse deflections of at least twice the plate thickness. As expected, the 
linear solution (which omits the higher order terms in the bracket) agrees with the 
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nonlinear solutions to transverse deflections of only about four-tenths of the plate 
thickness. 

§5.4 Change in Total Potential Energy of a Plate Without Delaminations 

Consider an isotropic, clamped circular aluminum plate of radius a subjected 
to a transverse load as in the previous section. The change in total potential energy 
II of the plate for an increment of the area (A = ira?) is 

G = -— — 

da 2ira 

where self-similar growth has been assumed. The change in total potential energy 
was calculated by DELAM as a function of the applied uniform load. For compar- 
ison, G was calculated using the perturbation method solution of Chia [41]. The 
toted potentied energy II was calculated using the stresses, strains, emd displace- 
ments as given by Chia (Appendix H), and differentiated with respect to the area 
A to give G. The present solution is compared to that of Chia as a function of the 
applied load in Figure 5-4. At lower loads, the solutions agree well. At higher loads, 
the perturbation solution for G is somewhat higher than the present solution, due 
to the use of more terms for the displacement functions in the present method. 
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Figure 5-3 Center deflection of a circular al uminum plate subjected to a uniform 
applied load (a = 1.0 in., h = 0.01 in.). 




Uniform Load, q^/Eh* 

Figure 5-4 Change in total potential energy of a clamped, circular aluminum 
plate subjected to a uniform load (a = 1.0 in., h = 0.01 in.). 
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Delaminations 

Several investigators have proposed approximate solutions to describe the be- 
havior of plates containing elliptical delaminations. To compare the present sublam- 
inate behavior model, which is quite general, to results presented in the literature, 
it is necessary to make several simplifications: 

1. The axes of the ellipse are aligned with the load axes. 

2. The sublaminate is a single layer which is either isotropic or orthotropic. 

3. The base plate is isotropic, and is much thicker than the sublaminate. 

The normalized critical buckling strain is defined as 

b 2 

€ « — £cr(l ~ V *V V V*)~77Ti \2 (®*®) 

\ h * ) 

where e cr is the far field strain in the plate when the sublaminate buckles, and h is 
the sublaminate thickness. The normalized critical buckling strain wits calculated 
by DELAM as a function of the plate ellipticity ( a/b ) for three cases: (a) an isotropic 
aluminum sublaminate and base plate, (b) a unidirectional sublaminate with the 
fibers aligned in the load direction (case A), and (c) a unidirectional sublaminate 
with the fibers aligned transversely to the load direction (case B). The base plate 
for these cases is a fictitious isotropic material. The material properties of the 
sublaminates and base plate are given in Table 5-1. 

The DELAM predictions axe compared with those of Chai and Babcock [18] 
and Kassapoglou [19] for the isotropic, orthotropic A, and orthotropic B sublami- 
nates in Figures 5-5 through 5-7, respectively. For the isotropic and orthotropic A 
sublaminates, the present results agree well with those of Chai and Babcock and 
Kassapoglou. For the orthotropic B sublaminate, the present results agree with 
Chai and Babcock and Kassapoglou at ellipticities greater than three. At lower 
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ellipticities, all of the analyses show different results. The orthotropic sublaminate 
B is very stiff compared to the base plate in the direction transverse to the applied 
load. Thus, under compressive loading the transverse Poisson expansion of the base 
plate drives the sublaminate into tension in the fiber direction. Conversely, an ap- 
plied tensile load will cause a Poisson contraction of the base plate and compression 
in the fiber direction of the sublaminate. Both the present solution and Chai and 
Babcock actually predict buckling under an applied tensile load (not shown) for 
case B at lower ellipticities. In any event, case B is an extreme situation for which 
none of the methods presently agree. 


Table 5-1 Material Properties used in the Comparisons 


Material 

Property 

Isotropic 

Aluminum 

Orthotropic 
Case A 

Orthotropic 
Case B 

Isotropic Base Plate 
for Cases A and B 

E x 

10.0 

1.47 

25.9 

1.47 

E y 

10.0 

25.9 

1.47 

1.47 

V XV 

0.30 

0.28 

0.016 

0.30 

G X y 

3.84 

1.03 

1.03 

0.567 

G xz 

3.84 

1.03 

0.286 

0.567 

Gy X 

3.84 

0.286 

1.03 

0.567 


The elastic moduli are in Msi. Poisson’s ratio is dimensionless. 
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§5.6 Summary 

The results indicate that the present analysis method can predict the buckling, 
postbuckling large deflection, and growth behavior of circular and elliptical plates 
and sublaminates subjected to various loads. Final verification of the method for 
the general cases of: (a) elliptical delaminations with axes arbitrarily oriented with 
respect to the applied in-plane loads, and (b) delaminations located between any 
two plies of the composite plate, will be made in Chapter 7 by comparison to 
experimental data. 
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Figure 5-6 Normalized critical buckling strain. Unidirectional type A sublam- 
inate on an isotropic base plate (Table 5-1) (6 = 0.5 in., h = 0.03 
in.). 



Figure 5-7 Normalized critical buckling strain. Unidirectional type B sublami- 
nate on an isotropic base plate (Table 5-1) ( b = 0.5 in., h = 0.03 
in.). 
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Experimental Procedure 


§6.1 Specimen Design and Fabrication 

A testing program provides experimented data against which the analytical 
development can easily be verified. The requirements for this experimental program 
were: (a) that the test specimens contain well-characterized delaminations; (b) that 
the specimens be exposed to uniform loads, implying a uniform far field strain; and 
(c) that the sublaminate deformation and growth be closely monitored. 

A sandwich construction test specimen (Figure 6-1) was developed of two 
Fiberite T300/976 graphite/epoxy face sheets secondarily bonded to an 0.625” thick 
aluminum honeycomb core. The honeycomb sandwich construction provided a test 
specimen that could be easily loaded in compression without introducing signif- 
icant bending moments. One of the face sheets contained a 0.001” thick Teflon 
disk between two plies simulating the presence of a delamination. The facesheets 
were fabricated from unidirectional T300/976 prepreg tape and cured in an auto- 
clave at a maximum temperature of 350° F at 80 psi. The secondary bonding of 
the facesheet laminates to the honeycomb was accomplished using a Hysol 250°F 
curing film adhesive under 30 psi in the autoclave. 
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The experimental parameters were delamination shape (circular or elliptical), 
orientation of the ellipse axes, delamination depth in the facesheet, and sublaminate 
layup. The specimens were nominally 3” wide by 6” long, and were fabricated 
in six groups, designated Series 1 through 6. After trimming, the ends of the 
specimens were filled with epoxy potting compound and milled flat and parallel to 
one another in preparation for testing. Series 1 through 3 were devoted to specimen 
development. Series 4 through 6 comprise the test matrix, the details of which are 
shown in Table 6-1. 

§6.2 Nondestructive Inspection 

Complete characterization of the delamination required an exact determination 
of the Teflon insert location in the facesheet. Despite careful positioning of the in- 
serts during fabrication of the laminates, trimming and milling operations changed 
the reference points. Every specimen was therefore ultrasonically C-scanned (Ap- 
pendix I) and the position of the Teflon insert mapped relative to the final dimen- 
sions. The location data were essential for the later mounting of the strain gauges. 
The C-scan dimensions of the delaminations were often 0.1” larger that the nominal 
size of the Teflon inserts, probably due to two phenomena: (a) incomplete bonding 
of adjacent plies at the edge of the Teflon, and (b) the lateral resolution of the 
C-scan. 

The C-scan was capable of mapping out not only the planar extent of the 
delamination, but also the depth of the delamination. This was useful as a check on 
the specimen fabrication, and particularly in the posttest inspection to determine 
whether delamination growth had occurred within the original ply interface or had 
progressed to other ply interfaces. 
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Table 6-1 Test Matrix 

Specimen 

Layup 

Teflon 

Orientation 

Depth* (plies) 

4-1 

[OieifOie] 

2” circle 

0° 

2 

4-2 

[0 16 if0 16 ] 

2” circle 

0° 

4 

4-3 

[OieifOie] 

2” circle 

0° 

6 

4-4 

[0i«J5T0 16 ] 

2” circle 

0° 

8 

5-1 

[(0290202902),if(sym)] 

2” circle 

0° 

3 

5-2 

[(0290202902)».ff(s!/m)] 

2” circle 

0° 

4 

5-3 

[(0290 2 02902),tf(sym)] 

2” circle 

0° 

5 

5-4 

[(0290202902)*/f(5ym)] 

2” circle 

0° 

8 

6-1 

[(O 29 O 2 O 29 O 2 )*iH r (5ym)] 

2” x 1.5” ellipse 

0° 

4 

6-2 

[(0290202902)#^T(aym)] 

2” x 1.5” ellipse 

30° 

4 

6-3 

[(O 2 902 O 29 O 2 ),H ( sym)] 

2” x 1.5” ellipse 

60° 

4 

6-4 

[(0290202902)*.ff(aym)] 

2” x 1.5” ellipse 

90° 

4 


* Number of plies from the facesheet surface. 


§6.3 Instrumentation 

Each specimen was instrumented with nine strain gauges arranged as shown 
in Figure 6-2. The individual type, orientation, and purpose of each gauge are 
shown in Table 6-2. In general, a single gauge was located at the center of the 
sublaminate to record the buckling and postbuckling strains during a test. One 
additional gauge was mounted on the front and three on the back facesheets to 
measure the far field strains in the specimen. The outputs of these gauges were 
used to adjust the relative load distribution between the faces during setup, to 
determine the actual load distribution during a test, and as a check against the 
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Figure 6-2 Specimen strain gauge locations and orientations. 

material properties given by the prepreg manufacturers . Four gauges were mounted 
on the periphery of the sublaminate, based on the C-scan data, to determine the 
onset of delamination growth. Up to the onset of growth, these gauges were also 
used to check the uniformity of the far field strains over the specimen. 



Sublaminate Side 
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Each specimen was loaded in compression between plattens in an MTS testing 
machine. One platten was a ball and socket self-aligning fixture to ensure that the 
loads were evenly distributed over the specimen. During each test, the plate was 
loaded at a constant displacement rate of .003 in/min. The outputs of all nine strain 
gauges (after amplification) and the MTS load cell were digitized and recorded in 
a spreadsheet computer file for later data reduction and plotting. Buckling of the 
sublaminate was observed, both visually and from the output of the strain gauge 
located at the center of the sublaminate. Growth of the delamination was detected 
by the four gauges surrounding the delamination. At extreme loads growth was also 
visually observed. 


Table 6-2 Strain Gauge Locations and Purposes 


Gauge Number 

Gauge Type 

Orientation* 

Purpose 

1 

CEA-06-062UW-350 

longitudinal 

sublaminate strain 

2 

CEA-06-125UN-350 

longitudinal 

fax field - front 

3 

CEA-06-125UN-350 

longitudinal 

fax field - back 

4 

CEA-06-125UN-350 

longitudinal 

far field - back 

5 

CEA-06- 1 25UN- 350 

longitudinal 

delamination growth 

6 

CEA-06- 1 25UN- 350 

transverse 

delamination growth 

7 

CEA-06- 125UN- 350 

longitudinal 

delamination growth 

8 

CEA-06- 125UN- 350 

transverse 

delamination growth 

9 

CEA-06- 1 25UN- 350 

transverse 

fax field - back 


* With respect to the load direction. 
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Comparison of Experimental and Model Results 


§7.1 Introduction 

This chapter demonstrates the validity of the delamination behavior model 
through a comparison with experimental data. Two kinds of data are required for 
this validation: (a) load-strain histories of delaminated sublaminates from the onset 
of loading through buckling and into the postbuckling regime, and (b) the load at 
which growth of the sublaminate begins. 

The experimental data in the literature can be divided into two types. In the 
first, experimental investigations have demonstrated the reduction in strength in 
composite plates resulting from impact damage. Data from a number of researchers 
have been reviewed by Baker et al. [42]. These data are important in that they 
were the first to show that a significant problem existed. However, the behavior 
of the sublaminate was not characterized in any way. In the second, researchers 
implanted a release agent, such as a teflon film, in the laminate during fabrication 
to simulate the presence of a delamination. Thus, a flaw of known shape, size, 
orientation, and depth in the laminate was introduced. Gillespie and Pipes [4], Wang 
et al. [6], Williams et al. [8], Wang and Slomiana [43], and Ramkumar [44] have 
simulated delaminations using through-width implants in wide columns. Whitcomb 


55 



Chapter 7: Comparison of Experimental and Model Results 


56 


[22], Wang and Slomiana [43], Ramkumar [44], Konishi and Johnston [45], Byers 
[46], Webster [47], and Geier et al. [48] have simulated delaminations using circular 
implants, while Jones et al. [17] and Reddy et al. [49] used rectangular implants and 
Kassapoglou [19] used elliptical implants. All have reported some features of the 
sublaminate behavior, such as buckling loads or growth loads. None, however, has 
reported a complete load-strain history. Therefore, an experimental investigation 
using implants w as undertaken to generate a data base for validation of the model. 

§7.2 Experimental Measurements and Material Properties 

Load-strain histories during compression testing were recorded from each of the 
nine strain gauges mounted on each specimen. The strain data from gauge number 
1, located at the center of the delamination, was used to establish the buckling 
and postbuckling behavior of the sublaminate. A typical response is illustrated in 
Figure 7-1. The onset of delamination growth and the corresponding growth load 
were determined from gauges 5, 6, 7, and 8, which were located on the periphery 
of the delamination. Typical responses are shown in Figure 7-2. As gauges 5 and 7 
began to deviate significantly from a linear response, the delamination had grown to 
include the gauges as part of the larger sublaminate. The growth load was estimated 
from the first gauge to show a change in linear behavior. 
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Figure 7*1 Typical sublaminate gauge 1 load-strain response. 



Figure 7-2 Typical sublaminate peripheral gauge 5, 6, 7, 8 load-strain responses. 
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The material properties of the cured laminate are shown in Table 7-1. These 
material properties were used in all analytical predictions. The compressive load 
was applied along the longitudinal xy axis of the specimens and the sublaminate 
orientations were measured counterclockwise with respect to this axis. The fol- 
lowing input data were common to all analyses (unless specified otherwise): (a) 
temperature difference, AT = — ISO 0 ! 71 ; (b) transverse pressure difference, A P = 0 
psi; and (c) relative growth direction, = 0. The strains shown in the comparisons 
were calculated at the outer surface of the sublaminate plus 0.004 in. to allow for 
the strain gauge thickness (xi» = 0., i2< = 0., xy — 4£- + 0.004). 


Table 7-1 Material Properties of Fiberite T300/976 


Material Property 

Value 

Units 

Longitudinal Young’s modulus, E x 

19.5E6 

psi 

Transverse Young’s modulus, E y 

1.32E6 

psi 

Poisson’s ratio, u xy 

0.30 

- 

In-plane shear modulus, G xy 

1.01E6 

psi 

Out-of-plane shear modulus, G xz 

1.01E6 

psi 

Out-of-plane shear modulus, G yz 

0.50E6 

psi 

Longitudinal thermal coeff. of expansion, a x 

0.50E-6 

in ^ 

Transverse thermal coeff. of expansion, a* 

18.0E-6 

in 

Critical strain energy release rate, G c [50] 

0.2 

in-lbf 

in 3 

Foundation modulus, K 

1.E6 

tn s 
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§7.3 Circular Delaminations in Unidirectional Laminated Plates 

Test Series 4 investigated the effect of delamination depth on the sublaminate 
behavior in unidirectional laminated plates containing circular delaminations. Tests 
4-1 through 4-4 used delamination depths of 2, 4, 6, and 8 plies, respectively. A 
small hole (0.021” diameter) was drilled through the sublaminate to allow air ingress 
to the Teflon implant to eliminate the effect of a transverse pressure differential. 
During the tests, specimens 4-1 through 4-3 were observed to buckle at increasing 
loads, while specimen 4-4 was loaded to the limit of the testing machine without 
buckling. Figure 7-3 shows the measured strains at the sublaminate center (gauge 
1) from each experiment. The predicted behavior of each experiment is shown as 
a solid line for comparison. The measured ply thickness (t = 0.00556”) was used 
as specific input for the analyses. Experimental and model results agree quite well. 
The data from the strain gauges surrounding the sublaminate (gauges 5, 6, 7, 8), 
indicating the onset of sublaminate growth, are shown in Figure 7-4. Only specimen 
4-2 showed an onset of sublaminate growth, which occurred at the end of the test. 

§7.4 Circular Delaminations in Cross Ply Laminated Plates 

Test Series 5 investigated the effect of delamination depth on sublaminate be- 
havior in cross ply ([(O 2 / 9 O 2 /O 2 / 9 O 2 )*# (sym)]) laminated plates containing circular 
delaminations. Tests 5-1 through 5-4 used delamination depths of 3, 4, 5, and 8 
plies, respectively. Specimens 5-1 ([02/90] sublaminate) and 5-3 ([O 2 / 9 O 2 /O] sub- 
laminate) were specifically intended to investigate the residual thermal strain effect, 
since the sublaminate layups were significantly different from the facesheet layup. 
A small hole (0.021” diameter) was again drilled through the sublaminate to allow 
air ingress to the Teflon implant. During the test, specimens 5-1 through 5-3 were 
observed to buckle at increasingly higher loads, while specimen 5-4 was loaded to 
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Figure 7-3 Load versus strain for a [0 ie /f0ie] T300/976 laminate under uniaxial compression. Circular de- 
laminations implanted at 2 (Specimen 4-1), 4 (Specimen 4-2), 6 (Specimen 4-3), and 8 (Specimen 
4-4) plies from the outer surface. Strain gauge 1 was located on the surface at the center of the 
sublaminate (see also Tables 6-1 and 6-2 and Figures 6-1 and 6-2). 
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Figure 7-4 Load versus strain for a [0i«/70i«] T300/976 laminate under uniaxial compression. Circular de- 
laminations implanted at 2 (Specimen 4-1), 4 (Specimen 4-2), 6 (Specimen 4-3), and 8 (Specimen 
4-4) plies from the outer surface. Strain gauges 5, 6, 7, and 8 were located around the periphery 
of the sublaminate (see also Tables 6-1 and 6-2 and Figures 6-1 and 6-2). 
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the limit of the testing machine without buckling. Figure 7-5 shows the measured 
strains at the sublaminate center (gauge 1) from each experiment. The predicted 
behavior of each experiment is shown as a solid line for comparison. The lami- 
nate cured ply thickness (t = 0.00609”) was measured and used as specific input 
for the analyses. The buckled region of each specimen in this series appeared not 
to extend over the full Teflon implant area but rather to have a shorter buckling 
dimension in the loading direction. The analyses confirmed that these specimens 
would preferentially buckle in multiple half waves in the loading direction and a 
single half wave in the tranverse direction. The analyses of specimens 5-1 and 5-3 
also indicated a small but noticable effect of the thermally induced load, in this 
case reducing the buckling load. The data from the strain gauges surrounding the 
sublaminate (gauges 5, 6, 7, 8), indicating the onset of sublaminate growth, are 
shown in Figure 7-6. Specimens 5-1 through 5-3 experienced sublaminate growth 
transverse to the applied load as indicated by peripheral gauges 5 and 7. 

§7.5 Elliptical Delaminations in Cross Ply Laminated Plates 

Test Series 6 investigated the effect of delamination orientation on sublaminate 
behavior for cross ply ([(02/90a/02/902)j7f(sym)]) laminated plates containing el- 
liptcial delaminations. Tests 6-1 through 6-4 used 2.0” by 1.5” elliptical Teflon 
implants oriented at 0°, 30°, 60°, and 90° to the applied load, respectively. Each 
implant was 4 plies deep in the facesheet. No hole was drilled through the sub- 
laminate so that any effect of transverse pressure might be observed. Buckling 
was observed in each specimen, followed by sublaminate growth. As in Series 5, 
the buckled region of each specimen appeared not to extend over the full Teflon 
implant area but rather to have a shorter buckling dimension in the loading direc- 
tion, which was confirmed by analysis. Figure 7-7 shows the measured strains at 
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Figure 7-6 Load versus strain for a [((h/902/(h/902)t If (0j/9(h/0W902)«] T300/976 laminate under uniaxial 
compression. Circular delaminations implanted at 3 (Specimen 5-1), 4 (Specimen 5-2), 5 (Spec- 
imen 5-3), and 8 (Specimen 5-4) plies from the outer surface. Strain gauges 5, 6, 7, and 8 were 
located around the periphery of the sublaminate (see also Tables 6-1 and 6-2 and Figures 6-1 
and 6-2). 







Chapter 7: Comparison of Experimental and Model Results 


65 


the sublaminate center (gauge 1) from each experiment. The predicted behavior of 
each experiment is shown as a solid line for comparison. The laminate cured ply 
thickness (t = 0.00609”) was measured and used as specific input for the analy- 
ses. Experiment and model results generally agree quite well. The data from the 
strain gauges surrounding the sublaminate (gauges 5, 6, 7, 8), indicating the onset 
of sublaminate growth, are shown in Figure 7-8. The observed growth direction 
was transverse to the applied load as indicated by peripheral gauges 5 and 7. The 
predicted growth load for specimens 6-1 and 6-2 was based on a relative growth 
parameter of = 0, while for specimens 6-3 and 6-4, it was based on = 100. 

§7.6 Buckling and Growth Loads 

The measured buckling loads are compared to the buckling loads predicted 
by the nonlinear equilibrium method in Figure 7-9. The measured and predicted 
growth loads are shown in Figure 7-10. The figures include the experimental and 
prediction uncertainties (Appendix J). The error bars shown are plus and minus 
three standard deviations. The dashed line in each figure represents perfect agree- 
ment between experiment and calculation. 

The buckling results show generally good agreement between prediction and 
experiment. The growth results show reasonable agreement. The prediction errors 
for the growth results are large due to large uncertainty in the critical strain energy 
release rate. Considering the uncertainties in the experimental data and the uncer- 
tainties in the analyses, it would appear that the analysis method describes with 
reasonable accuracy the experimental data. 
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Figure 7-8 Load versus strain for a [(0 2 /90 2 /(h /90 2 L //(0 2 / 9 O 2 /O 2 /90 2 ) , ] T300/976 laminate under uniaxial 
compression. 2.0 in. by 1.5 in. elliptical delaminations implanted 4 plies from the outer surface 
(all specimens). 0° (Specimen 6-1), 30° (Specimen 6-2), 60° (Specimen 6-3), and 90° (Specimen 
6-4) orientations of the ellipses with respect to the load. Strain gauges 5, 6, 7, and 8 were located 
around the periphery of the sublaminate (see also Tables 6-1 and 6-2 and Figures 6-1 and 6-2). 
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re 7-9 Measured versus predicted buckling loads. Error bars are plus 
minus three standard deviations (see Tables J-l and J-3). 



Predicted growth load, N g Qbf/m) 

Figure 7-10 Measured versus predicted growth loads. Error bars are plus and 
minus three standard deviations (see Tables J-l and J-3). 






Chapter 8 


Sample Problem and Discussion 


§8.1 Introduction 

The analytical and experimental verifications performed in Chapters 5 through 
7 demonstrated the model performance over a range of input variables and load- 
ing conditions. The uncertainty analysis (Appendix J) determined not only the 
overall uncertainty in the experimental data and code predictions, but the relative 
sensitivity of the model to specific input variables. In general, the input variables 
can be grouped into three types in terms of their influence on the buckling load, 
postbuckling strain, and growth load (Table 8-1): (a) variables which are physically 
well-characterized and for which the results are highly sensitive, (b) variables which 
are physically well-characterized but for which the results show little sensitivity, 
and (c) variables which, for a variety of reasons, are poorly characterized and thus 
may exert a large influence on the results. This chapter discusses the effects of this 
last group of variables in the context of a sample problem. 

§8.2 Sample Problem Description 

The sample problem chosen corresponds to Experiment 6-1. The plate is a 
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16 ply symmetric cross ply laminate ([02/902/02/902],) fabricated from Fiberite 
T300/976 graphite/epoxy (see Table 7-1 for material properties). A 2.0 in. by 1.5 
in. elliptical delamination lies 4 plies deep in the plate. The ellipse major axis is 
aligned with the load axis. The plate is subjected to uniform compression along the 
major axis of the ellipse as shown in Figure 8-1. A complete list of the nominal input 
variables is shown in Table 8-2. Note that xj, X 2 , and X 3 are the coordinates in the 
sublaminate, and in this instance are coincident with the xi», xy, x$> coordinate 
system of the plate (see Figure 3-1). The strains shown in the examples were 
calculated at the outer surface of the sublaminate plus 0.004 in. to allow for the 
strain gauge thickness (xi = 0., X 2 = 0., X 3 = Ap + 0.004). 


Table 8-1 Code Input Variable Sensitivities 



Sensitive 

Insensitive 

Well- 

Ellipse major semi-axis, a 

Long. Young’s modulus, E x 

Characterized 

Ellipse minor semi-axis, b 

Tran. Young’s modulus, E y 


Ellipse orientation, Q 

Shear modulus, G xy 


Ply thickness, U 

Long, thermal expansion, a x 
Tran, thermal expansion, a y 
Temperature change, AT 

Poorly 

Critical strain energy release, G c 


Characterized 

Growth direction parameter, || 
Foundation modulus, K 



Transverse pressure load, A P 
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In each example that follows, one variable has been allowed to change over a 
realistic range of possible values to illustrate the effect of that variable while holding 
all other input variables constant. The effects of each variable on the calculated 
response are shown individually in the following sections. Recommendations for 
designers are summarized in the last section. 
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Figure 8-1 Sample problem description. Input variables sire given in Table 8-2. 
Material properties are given in Table 7-1. 


§8.3 Geometry Effects 

It is well known from classical buckling theory that the buckling load of the 
sublaminate varies in proportion to the cube of the thickness and inversely with 
the square of the lateral dimensions of the sublaminate. For a designer studying 
the effects of manufacturing-induced disbonds or impact-caused delaminations, the 
thickness and shape of the sublaminate may be the source of major uncertainty. 
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The effect of changing the nominal ply thickness in the sample problem on the 
calculated load-strain history of the sublaminate is illustrated in Figure 8 - 2 . The 
three cases represent successive changes of 0.0002 inch to the ply thickness. The 
buckling and growth loads are shown versus the ply thickness in Figure 8 - 3 . The 
buckling load and growth load increase by roughly 100 lbf/in with the total 0.0005 
in. thickness variation. 


Table 8-2 Sample Problem Input Variables 


Variable 

Value 

Units 

Material 

T300/976 


Layup 

(0290202902)*Ff(5ym) 


Delamination depth 

4 

plies 

Ellipse major semi-axis, a 

1.00 

in 

Ellipse minor semi-axis, b 

0.75 

in 

Ellipse orientation, 9 

0. 

degrees 

Ply thickness, t, 

0.00556 

in 

Critical strain energy release rate, G c 

0.20 

tn T 

Relative growth direction parameter, jffc 

0 


Contact law foundation modulus, K 

1 .E 6 

4*4 

»n 3 

Normal load in the 1 direction, 71 

1 . 


Normal load in the 2 direction, 72 

0 . 


Shear load in the 1-2 plane, 7 6 

0. 


Change from reference temperature state, AT 

-180. 

°F 

Transverse pressure load, A P 

3. 

psi 




Applied Load, N (lbf/in) «g‘ Applied Load, N (lbtfin) 
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2000 


ItXIl 


Strain, e i (10*” in/in) 


ire 8-2 Effect of changing ply thickness on the compressive load-strain re- 
sponse of the sample problem described in Figure 8-1. Results calcu- 
lated by the DELAM code. 


3000 


2000 


Ply Thickness, t (m) 

Figure 8-3 Effect of changing ply thickness on the buckling and growth loads of 
the sample problem described in Figure 8-1. Results calculated by 
the DELAM code. 
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The effect of a change in the lateral dimension of the sublaminate on the load- 
strain behavior is more dramatic than the effect of a thickness change. The effect 
of 0.05 inch successive changes in the transverse “b” dimension of the delamination 
ellipse on the load-strain behavior of the sample problem sublaminate is shown in 
Figure 8-4. The larger diameter sublaminates are much more compliant than the 
smaller sublaminates. The trends of the buckling and growth loads versus the ellipse 
semi-axis dimension sire shown in Figure 8-5. As expected, the loads decrease with 
an increase in the semi-axis. These effects are particularly important because in 
practice, the lateral dimensions of a delamination may only be known to about the 
accuracy shown in this figure. 

§8.4 Contact Model Effects 

The contact model represents the physical restraint to deflection of the sub- 
laminate posed by the plate containing the delamination. The key to the model 
is the value of the foundation modulus K (Equation 3.31). No value for K has 
been measured for graphite/epoxy. A rationale for estimating K from the t ran verse 
elastic modulus of the foundation Ef and a characteristic length If is discussed in 
Appendix D. For the materials used here, Ef « 1.0 Msi and If « 1.0 inch yielding 
a foundation modulus K « l.xlO -6 (see Equation 3.32). Figure 8-6 illustrates 
the effect on the load strain history of the sublaminate of successive changes in the 
foundation modulus K. The value of K = 0 indicates that the contact model was 
not used. An incease in the foundation modulus corresonds to an increase in the 
stiffness of the response. Varying the value of the foundation modulus has almost no 
effect on the buckling and growth loads of the sublaminate, as shown in Figure 8-7. 



Applied Load, N (lbf/in) <w* Applied Load, N (lbf/in) 
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3000 


b = 0.65 in 


2000 


Strain, £j (10’° in/in) 

ire 8-4 Effect of changing the semi-minor axis “b” of the ellipse on the com 
pressive load-strain response of the sample problem described in Fig 
ure 8-1. Results calculated by the DELAM code. 


3000 


2000 


Ellipse Semi-Axis, b (m) 

Figure 8-5 Effect of changing the semi-minor axis “b” of the ellipse on the buck 
ling and growth loads of the sample problem described in Figure 8-1 
Results calculated by the DELAM code. 
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3000 r K = 10 7 lbf/in 3 


2000 


Strain, e; 


ln/in) 


Figure 8-0 Effect of changing the contact law foundation modulus K on the com 
pressive load-strain response of the sample problem described in Fig 
ure 8-1. Results calculated by the DELAM code. 


Foundation Modulus, K (lbf/in-’) 

Figure 8-7 Effect of changing the contact law foundation modulus K on the buck- 
ling and growth loads of the sample problem described in Figure 8-1. 
Results calculated by the DELAM code. 
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The transverse pressure model describes the effects of subatmospheric pres- 
sure in the cavity formed by the sublaminate as it buckles away from the plate. 
Since there is no method to measure the actual A P across the sublaminate, the 
uncertainty associated with A P may be large. Figure 8-8 illustrates the effect of 
the pressure differential on the load-strain response of the sublaminate. Figure 8-9 
shows the increasing buckling and growth loads with increasing pressure differential. 

§8.6 Growth Model Effects 

The growth model requires a parameter, describing the shape of sublami- 
nate growth, and a material property, G c , which is the critical strain energy release 
rate of the material. Neither is well-characterized [50]. The effect of changing values 
of on the calculated growth load of the sample problem is illustrated in Figure 8- 
10. The lowest value is clearly ^y = 0, and the growth load increases sharply with 
increasing values of ^y. Figure 8-11 depicts the dependence of the growth load on 
the critical strain energy release rate. As shown, the growth load increases strongly 
with increasing G c < 
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Figure 8-8 Effect of changing the transverse pressure A P on the compressive 
load-strain response of the sample problem described in Figure 8-1. 
Results calculated by the DELAM code. 



Figure 8-9 Effect of changing the transverse pressure A P on the buckling and 
growth loads of the sample problem described in Figure 8-1. Results 
calculated by the DELAM code. 
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Figure 8-10 Effect of changing the growth model parameter on the compressive 

load-strain response of the sample problem described in Figure 8-1. 
Results calculated by the DELAM code. 



Critical Strain Energy Release Rate, G c flbf-in/in^) 

Figure 8-11 Effect of changing the critical strain energy release rate G c on the 
growth loads cm the sample problem described in Figure 8-1. Results 
calculated by the DELAM code. 
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§8.7 Summary and Recommendations 

The sample problem shown has demonstrated the sensitivity of the model to 
certain input variables, in particular the sublaminate geometry, the foundation mod- 
ulus, the transverse pressure differential, and the growth parameter and critical 
strain energy release rate. In summary, none of the variables affect the prebuckling 
behavior of the sublaminate whereas all of them affect the postbuckling load-strain 
behavior. Of the variables affecting the buckling load and growth load, the sub- 
laminate thickness t and the foundation modulus K were shown to have a minor 
effect. In contrast, the sublaminate lateral dimension “b,” the pressure difference 
A P, and the growth parameter || and the critical strain energy release rate G c 
were all shown to significantly affect the buckling and growth loads. Therefore, 
the following choices are recommended to designers for conservative analysis: (a) 
the lateral dimensions should be chosen large, (b) the pressure difference should be 
zero (which is likely since the pressure can equalize through cracks), (c) the growth 
parameter ^ should be systematically evaluated to find the lowest growth load, 
and (d) the critical strain energy release rate should be as low as practical. 
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Concluding Remarks 


A model was developed to describe the behavior of delaminated composite 
plates subjected to compressive in-plane loads. The delaminated region is assumed 
to be elliptical, and may be located between any two plies of the laminate. The axes 
of the ellipse may be arbitrarily oriented with respect to the applied loads. The 
model calculates the displacements, strains, and stresses in the plate containing 
the delamination, and in the sublaminate created by the delamination. The model 
solves the nonlinear equilibrium equations describing the sublaminate up through 
large postbuckling deflections of the sublaminate. In particular, the model pre- 
dicts the loads applied to the plate at which first buckling and then growth of the 
sublaminate will occur. 

A computationally efficient computer implementation of the model was devel- 
oped. The code has a user friendly interface, and is intended to be used for design 
calculations. 

A new set of experimental data on the behavior of Fiberite T300/976 graphite/ 
epoxy laminated plates containing simulated delaminations and loaded in compres- 
sion was used to validate the model performance. The sublaminate load-strain 
histories were described at a level of detail not previously available in the literature, 
and will prove useful in future delamination studies. 
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The model currently describes a single delamination, and the behavior of the 
sublaminate is assumed not to affect the behavior of the plate in which it is con- 
tained. The effects of multiple delaminations, and the interaction of the sublaminate 
and plate, are suitable topics for future investigations. 
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Appendix A 


Engineering Constants for Isotropic and Orthotropic Materials 


This appendix identifies the ply plane stress reduced stiffnesses Qij used in 
Chapter 3 in terms of ordinary engineering elastic constants for isotropic and or- 
thotropic materials. In addition, the transformation matrix relating on-axis and 
off-axis stiffnesses is given. 


§A.l Engineering Constants 


The on-axis lamina constitutive relations for an orthotropic material in plane 
stress are 


( cr *\ 


{ Qxx Qxy 0 0 0 ^ 


(€ X \ 

<7y 


Qxy Qyy 0 0 0 


€ v 
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0 0 Qpp 0 0 


e p 
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0 0 0 Qrr 0 
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where x and y are the in-plane principal material axes of the lamina, p refers to the 
transverse y-z plane, r to the transverse x-z plane, and s to the in-plane x-y plane. 


In terms of engineering constants, the stiffnesses for an orthotropic material axe 
Ex _ Ey 


Qxx — 


1 - VxyVyx 


Q P p — G p 


Qyy ~ 
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rr — G r 
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Qaa — G, 


1 - VxyVyx 


(A.2) 


where 
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E x = Longitudinal Young’s modulus 

E y = Transverse Young’s modulus 

u xy = Major Poisson’s ratio 

v yx = Minor Poisson’s ratio 

G p = Transverse shear modulus in the y-z plane 

G r = Transverse shear modulus in the x-z plane 

G a = Transverse shear modulus in the x-y plane 

For an isotropic material the relationships are simpler. In terms of engineering 
constants the stiffnesses for an isotropic material are 


Qxx = 2 QyV = 1 - 1/2 Qxy = 1^2 
Qpp = G Qrr = G Qgg = G 

where 

E = Young’s modulus 
v = Poisson’s ratio 
G = Shear modulus = 2 ( 1 + ») 


( 24 . 3 ) 


§A.2 Transformation Matrix 

The off-axis l amin a constitutive relations for an orthotropic material in plane 
stress are 


/<71\ 


/Qu Q 12 0 0 Qie\ 


(ti\ 
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Ql2 Q 22 0 0 C?26 
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0 0 Q 44 Q 45 0 


C4 
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0 0 Qm Qh 0 
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\<7 J 


\C?i6 Q 26 0 0 Qm/ 


\ej 


where 1 and 2 are the in-plane body axes of the lamina, 4 refers to the transverse 
2-3 plane, 5 to the transverse 1-3 plane, and 6 to the in-plane 1-2 plane. The 
transformation matrix relating the off-axis stiffnesses to the on-axis stiffnesses is 
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where m = cos8, n = sm0, and ^ is the angle between the on-axis and off-axis coor- 


dinate systems, defined as positive in the counter-clockwise direction. Equation A. 5 
is given for the negative transformation, meaning that the ply on-axis stiffnesses are 


rotated to the body off-axis coordinates. 


Appendix B 


Integration of the Plate Strain Expressions 


This appendix details the integration of the plate strain expressions (Equation 
3.12) given in Chapter 3. The strain field in the plate is constant and uniform over 
the plate. The strains were derived in terms of constants c, describing the plate 
(Equation 3.13), the load N, and the residual thermal strains <, e? lT (Equation 3.14): 



N + 



(B . 1 ) 


In terms of the displacements the strains are 



‘ \ 

ay 1 

dx\ 

\ 9 y*‘ j . / 

' SX2 dXi ' 


(B. 2) 


Combining Equations B.l and B.2 and integrating the first two expressions yields 


°< = (ciiV + + /(* 2) 

%i' = (c 2 N + ’’4 IT )x 1 +g(t 1 ) 


(B.3) 


where / and g are arbitrary functions of x% and xj, respectively. The in-plane shear 
strain °^ 1 can then be expressed as 


92 



Appendix B: Integration of the Plate Strain Expressions 


93 


, _ ax' . ax' _ 4F(* 2 ) . M*i) _ r v . v it 

e 6 - “31“ + IT" - ~JZ — + —J1 — - c « iV + € e 


(5.4) 


dxi ' dx\ dx 2 dxi 

Since the strains in Equation B.l are constant, f(x 2 ) and gr(xi) can be at most 
linear functions of x 2 and Xj, respectively. Assuming that no rigid body rotation of 
the plate occurs, then [54] 

(B.5) 

dx 2 dx\ 

or, combining Equations B.4 and B.5 




(5.6) 


Finally, and are easily integrated and substituted into Equation B.3 to 

yield the in-plane displacements in the plate 


M 1 = + °^ ,t )®2 

°u$' = ^(cgJV + ° e g ,T )r 1 + (c 2 JV + °e5 ,T )x 2 


(5.7) 


Appendix C 


Basic Assumptions of Nonlinear Plate Theory 


This appendix states the basic assumptions of the nonlinear plate theory used 
in the main text. In particular, the plausibility of some of the assumptions is 
demonstrated via an order of magnitude calculation. The motivation for this chapter 
was derived from the realization that the Kirchhoff-Love assumption (that normals 
to the midsurface remain normal) was an integral assumption in von Karman’s 
nonlinear large deflection plate theory. The appropriate and consistent nonlinear 
strain measures for the shear deformable plate theory are developed subsequently 
[ 41 , 51 ], 

§C.l Basic Assumptions 

Given: A plate geometry. A body B is defined by two parallel surfaces and an 
edge surface joining them such that characteristic in-plane (xi and x 2 coordinates) 
lengths ( L ) are much greater than the through- the- thickness (£3) length (h). 


h<cL 


(C. 1) 


Assumption 1: Tangential displacements ui and u 2 are infinitesimal but the trans- 
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verse displacement U3 is on the order of the plate thickness. 

ui, U2 ~ small 
U3 ^ 0 (h) 

Assumption 2: Derivatives of displacements are moderate to small. 

|i~0(£)orO(£)’ U= 1,2,3 

Assumption 3: Linear strain parameters e tJ are small, and linear rotation 
ters u>ij are moderate. 


e ij 


1 ,dui dui. . . 

iS&Tj + *,; = 1,2,3 


1 ,dui duj h. 

UiJ 2 ^dxj dx t,J ~ 1,2,3 


Assumption 4: In-plane rotations are small. 


U > 12 ~ 0(-) 2 


§C.2 Order of Magnitude Estimates 
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(C. 2) 


(C.3) 

parame- 

(C.4.a) 


(C.5.6) 


(C.6) 


Figure C-l Order of Magnitude Estimates for Displacements and Derivatives 
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To see why the above assumptions sire reasonable, consider a plate with length 
dimensions of L and thickness h, with a transverse displacement on the order of h 
at one end (Figure C-l). In the deformed configuration the plate will have been 
stretched, rotated, and thinned by this displacement with respect to the original 
configuration. With respect to a Cartesian coordinate system fixed in the original 
reference configuration (Lagrangian description), the change in tangential in- plane 
displacement with respect to the in-plane coordinate is of the order 

V TT+h?-L 

L 

Wi + Cx) 2 - 

~0(-* 

1j 

~0( l + i(£)* + . 

~<x I > 2 

Similarly, the change in transverse displacement with respect to the in-plane coor- 
dinate is of the order 


dui 

dxi 


0 ( 


L 


(C.7) 


.- 1 ) 


< c - 8 > 

The change in the transverse displacement with respect to the transverse coordinate 
due to rotation (the shear effect is smaller) is of the order 

*! i „ 

dx 3 v h } 

(a9) 


Similar estimates for the remaining derivatives can be made so that, in summary, 
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the derivative order of magnitudes are 
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( 0 . 10 ) 


as stated in (C.3). 

Using the estimates for the magnitudes of derivatives (C.10), the plausibility 
of Assumption 3 for linear strain and and rotation parameters can be checked. For 
example, 


0“i h 2 

e " = dT t ~ 0( 3> 


du 2 h 2 

” W, ~ 0{ l } 


du 3 h j 

= dT, ~ 0{ l) 


(C.ll) 


The linear shear strain parameters and the rotation parameters must be examined 
carefully. In particular, the algebraic values of the derivatives 8X6 opposite in 
sign to the algebraic values of for a rigid body rotation. Thus, the linear shear 
strain parameters are really difference equations and the rotations simply additive. 
Therefore, Assumption 3 actually states that the difference in derivatives is small. 
That is, 
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In the above, we implicitly used a Lagrangian description of the system. That 
is, within a fixed Cartesian coordinated system, the deformation of the plate was 
described in terms of the reference undeformed configuration of the plate. The 
appropriate stress measure in the Lagrangian description [52] is the 2nd Piola- 
Kirchhoff stress tensor which, in terms of the Cauchy stresses <r-, is 


r P-K _ Po dXj c 


(C M3) 


« P'dxidzf** 
where p 0 and p c axe the mass densities in the reference and current configurations, 
respectively, and x° and x% particle locations in the reference and current configu- 
rations. The appropriate strain measure [52] for the Lagrangian description is the 
Green-Lagrange strain tensor 


1 , duj duj_ du k duic . 

e,JI 2 dij dii dii dxj 


(C.14) 


The 2nd Piola-Kirchhoff stress tensor and the Green-Lagrange strain tensor 
axe energy conjugate to one another [53]. That is, the strain energy of the system 
calculated by these two measures in the reference coordinates is equal to the strain 
energy calculated in the current coordinates using Cauchy stresses and infinitesimal 
strains. Thus, we have defined a plate theory involving large transverse deflec- 
tions, moderate rotations, and small strains. For nonlinear analysis, the following 
observations can be made. 


C.3.1 Stresses 

The Cauchy stresses, expressed in terms of the 2nd Piola-Kirchhoff stresses and 
the derivatives of displacement already discussed, axe 

T C Pc r P-K , (C & u i | c & u j , & u j\ P-K} /'Cl'i'l 

'o-rK +{S ^^r + SiQ dr 0 + d^d^ )<r ^ } (C7 * 15) 
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where Sy is the Kronecker delta and tx, and x, are with respect to the reference 
configuration. Given the magnitudes of the displacement derivatives (C.10), the 
consequence of these assumptions is that, to first order, the Cauchy stresses and 
2nd P iola- Kir chhoff stresses are equal. Therefore, the stress and moment resultants 
defined elsewhere in this thesis will also be approximately equal. There is no further 
need to distinguish reference and current configurations when discussing stress. 

C.3.2 Strains 

The components of the Green-Lagrange strain tensor (C.14) and the order of 
magnitude of the various displacement derivatives are: 


, _ - , * /•/''-I \2 , I \2\ 

— + + (3^) +(— )) 


5x x 2"dxi 


dxi 


‘n~0{(±)\ (£)‘, (*) 4 , (h.f) 


(C. 16) 
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_ 1 du 2 duz dui dui du 2 du 2 duz du 3 

623 2 dxz dx 2 dx 2 dxz dx 2 dxz + dx 2 dxz 


<23 ~ 0 {( T ), ( 7 ), (~) 3 , ( T ) 3 , 


(jf) 


(C. 21) 


v £' ’ v £' 

Neglecting terms of (£) 3 and higher, and recalling the arguments for the magnitudes 
of the linear portions of the strains (C.12), it is apparent that 


1. ) Only the in-plane strains retain nonlinear terms, and 

2. ) With the exception of e 33 , these terms are derivatives of the transverse dis- 

placement Uz. 


In this regard e 33 deserves a special note. Recalling the estimate of the magni- 
tude of due to rotation (C.9), it is apparent that and are algebraically 
opposite (similar to the linear shear strain terms). This makes sense in that a strict 
rigid body rotation should produce no strain. For moderate rotation deformation 
the rigid body portion of that motion should vanish. The difference between these 
terms will be due to any effect of shear on e 33 (arguably of order (£) 3 ) and the 
term will be due to Poisson thinning. (These are obviously beam-like simplifi- 
cations.) Thus, a reasonable summation of appropriate nonlinear strain measures 
for the above assumptions in a compact form is 


_ 1 dm duj duzduz . .. 10 o 

Eii “ 2^dxj + dxi + dxi dx^ ,3 l * ’ 3 


(C.22) 


That is, the strain-displacement relations von Karman assumed are entirely appro- 
priate in the context of moderate rotation, small strain shear deformation theory. 
For large strains, particularly large shear strains, appropriate nonlinear terms will 
have to include the effects of the assumed shear deformation mode in the nonlin- 


earities. 




Appendix D 


Contact Model Foundation Modulus 


The contact model requires the foundation modulus K of the plate as a material 
property input. For graphite/epoxy composite materials no measured value of K 
is available. The purpose of this appendix is to describe a method for estimating 
this property. The vertical displacement v of a semi-infinite plate subjected to a 
uniform pressure load q acting on a portion of the plate, as shown in Figure D-l, is 
[54] 


2q 

v — -^2 a log a (.D.l) 

irE 

where the displacement is evaluated at the origin, E is Young’s modulus, and 2a is 
the width over which the pressure acts. Rearranging Equation D.l gives 


9 = 


kE 

4 a log a 


The contact model is stated as 


(D.2) 


, (A P- Ku\ l u\ l > 0 

/ = |a p 4'<o 


(D.3) 


where / is the contact force per unit area and u sl is the displacement of the sub- 
laminate. Identifying q with / and v with u sl , the foundation modulus K may be 
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estimated from 


K = 


wE 

4 a log a 


More simply, the foundation modulus is estimated by 


(DA) 



(D- 5) 


where Ej is the elastic modulus of the foundation and If a characteristic length. 



x 3 


Figure D-l Uniform pressure load acting on a semi- infini te plate. 


Appendix E 


Parallel Axes Theorem for Unsymmetric Laminates 


This appendix describes the calculation of laminate stiffnesses for generally 
unsymmetric laminates. Aij, Bij , Dij, Eij, Fij, Hij are laminate stiffnesses defined 
in Equation 3.34 as 


*r 

(Aij , Bij , Dij , Eij ,Fij, Hij ) = J Qij(l,x3,xl,xl,xl,xl)dx 3 i,j = 1,2,6 


V 

(Aij, Dij, Fij) = J Qij(l,xl,x^)dx3 i,j — 4,5 


(E.l) 


For a laminate composed of plies of varying thicknesses or an odd number of plies 
the laminate midsurface may fall within a given ply. It is computationally simple 
to calculate the laminate stiffnesses in a coordinate system originating on the lam- 
inate outer surface, and then to use the parallel axis theorem [23] to determine the 
stiffnesses in a coordinate system located at the laminate midsurface. 


For a primed coordinate system located a distance d from the laminate mid- 
surface as shown in Figure E-l, the laminate stiffnesses are calculated from 
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<-*r 

d+*£ 

= j Qij(\,x\,x*)dx' 3 i,j = 4,5 
d-*i 1 


= 1 , 2,6 


(E. 2) 


Notice that the stiffness G\j is required. The laminate stiffnesses in the imprimed 
plate midsurface coordinate system are then determined using the parallel axis 


theorem as 


A-- = A' 

Bij — Bjj dAij 

Dij « - 2 dBij - cPAij 

Eij = E' tj - ZdDij - 3d 2 Bij -d 3 Aij (E.3) 

Fa = F[j - 4 dEij - 6d?Dij - 4d 3 B tj - dM 0 

Gij = G'ij - 5 dFij - 10 #Eij - 10d 3 I> o - 5d*Bij - d 5 ^,, 

Hij = ^ - 6dGij - lSd 2 ^ - 20d 3 Eij - - 6d 5 B 0 - d 6 ^ 


Each equation requires the result of the previous equation to complete the calcula- 
tion. The distance d is, in general, arbitrary. However, in the particular case here, 
d is equal to half the sublaminate thickness . 
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Figure E-l Laminate thickness direction coordinate systems located at the plate 
midsurface (£3) and an arbitrary distance d from the plate midsurface 

(*i). 



Appendix F 


Strain Energy Release of an Elliptical Sublaminate 


This appendix details the derivation of the strain energy release of an elliptical 
sublaminate (Equations 3.41 and 3.43). The total potential energy II and area A 
of an elliptical plate are functions of the semi-major and semi-minor axes a and 6, 
respectively, of the ellipse. 

II = II(a, b , ...) A = irab (Fl) 


Taking differentials of both yields 


<m = da + ^db dA = ir(bda + ada ) 

Combining the differentials gives the strain energy release per unit area 

_ jgda + jfrdb 
dA ir(bda + ada) 


(F.2) 


(F.3) 


or 


jtt on da * on 

«5 = hjlIm 
+ a ) 


(F.4) 


In the particular instance of a sublaminate in a composite plate, the calculation is 
made for each system, IP* and IP*. 
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DELAM SAMPLE INPUT/OUTPUT 


THIS IS PROGRAM OELRM. 

COPVRIGHT 1989 BV SCOTT OUEN PECK 

STRUCTURES AND COMPOSITES LRBORRTORV 
DEPARTMENT OF AERONAUT I CS RNO ASTRONAUTICS 
STANFORD UNIUERSITY, STANFORD, CALIFORNIA 
(415) 723-4135 

GIUEN A LAMINATED COMPOSITE PLATE CONTAINING AN 
ELL I PT I CALLV SHAPED DELAMINATION, OELAM UILL 
CALCULATE THE FOLLOUING: 

(1) THE CRITICAL LOAD APPLIED TO THE PLATE NECESSARY 

TO CAUSE BUCKLING OF THE SUBLAMINATE CREATED BV 

THE DELPHI NAT I ON, 

(2) THE NONLINEAR LOAD-STRAIN HISTORY OF THE 

SUBLRMINRTE, RND 

(3) THE CRITICAL LOAD APPLIED TO THE PLATE NECESSARY 

TO CAUSE THE ONSET OF OELAMINATION GROUTH. 

I********************************************************************* 


DELAM IS BASED ON A NONLINEAR PLATE THEORY 
INCLUDING THE EFFECTS OF LARGE TRANSUERSE DEFLECTIONS 
OF THE SUBLAMINATE AND TRANSUERSE SHEAR DEFORMATION. 
THE ASSUMPTIONS ABOUT THE PLATE AND SUBLAM I NATE ARE: 


(1) THE PLATE CONTAINING THE DELAMINATION IS 
SYMMETRICALLY LAMINATED. 


(2) THE DELAMINATION MAY OCCUR BETUEEN ANY TUO 

PLIES, AND THEREFORE THE SUBLRMINRTE MAY 
BE UNSVHMETRICALLY LAMINATED. 

(3) THE ELLIPTICAL SUBLAMINATE MAY BE ARBITRARILY 
ORIENTED UITH RESPECT TO THE APPLIED LOADS. 

(4) THE PLATE FORMS A DETACHED ELASTIC FOUNDATION 

FOR THE SUBLAM INATE, UHICH MODELS POSSIBLE 
CONTACT BETUEEN THE TUO. 


(5) A TRANSUERSE PRESSURE DIFFERENTIAL HAY ACT 
ACROSS THE SUBLAMINATE THICKNESS DUE TO 
SUBATMOSPHERIC PRESSURES IN THE CRUITY FORMED 
BETUEEN THE SUBLRMINRTE AND PLATE. 

(6) GROUTH OF THE SUBLAMINATE UILL OCCUR 

UHEN THE TOTAL POTENTIAL ENERGY RELEASED BV 
SUBLRMINRTE - PLATE SYSTEM EXCEEDS THE 


CRITICAL UALUE FOR THE PARTICULAR MATERIAL. 
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DELAH OPERATES IN ONE OF TUO NOOES: PROflPTED 
AND ORTA FILE. IN THE PROflPTED MODE VOU UILL BE QUERIED 
FOR EACH INPUT PARANETER. IN THE ORTA FILE NODE, THE 
INPUT DATA IS ASSURED TO BE IN A USER OATA FILE, 

AND VOU UILL BE ASKEO ONLV FOR THE NANE OF THE FILE. 

AT THE ENO OF THE INPUT PROCESS IN PRONPTEO NOOE, 

VOU UILL BE ASKED UHETHER VOU UOULD LIKE THE INPUT 
TO BE SAUEO IN AN INPUT FILE FOR FUTURE ANALVSES. 

aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa 


UOULD VOU LIKE THE NORNAL PRINTOUT (0), 

MORE PRINTOUT (1), OR LOTS OF PRINTOUT (2) ? 

0 

UOULD VOU PREFER PRONPTEO INPUT (P) 

OR TO READ VOUR INPUT FRON A DATA FILE (D)? 

P 

aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa 


THE FOLLOUING INPUT DESCRIBES THE GEONETRV AND 
HATER I ALS OF THE ELLIPTICAL DELAN I NAT I ON . 

THE INPUT UNITS ARE IN ANV SELF-CONSISTENT 
SVSTEH THE USER OESIRES. 


ELLIPSE SEN I -MAJOR AXIS? 

1 .0 

ELLIPSE SEN I -H I NOR AXIS? 

0.75 

ROTATION OF THE ELLIPSE U.R.T. THE PLATE? 

0 . 

NUMBER OF PLIES IN THE UHOLE PLATE? 

1 6 

NUMBER OF PLIES IN THE SUBLAMINATE? 

1 

THE PLIES ARE NUMBERED FROM THE TOP SURFACE TO THE BOTTOM SURFACE. 
SHOULD EACH PLV HAUE THE SAME THICKNESS? (V/N) 

V 


PLV THICKNESS - 
.00556 

THE ORIENTATION OF EACH PLV IS POSITIUE FROM THE PLATE COORDINATE AXIS 
TO THE PLV AXIS. 

PLV NUMBER f ORIENTATION - 

0 . 

PLV NUMBER 2 ORIENTATION - 

0 . 




PLV 
90 , 

HUNGER 

3 

ORIENTATION 

PLV 

90. 

NUNBER 

4 

ORIENTATION 

PLV 
0 . 

NUNBER 

5 

ORIENTATION 

PLV 
0 . 

NUNBER 

6 

ORIENTATION 

PLV 

90. 

NUNBER 

7 

ORIENTATION 

PLV 

90. 

NUNBER 

9 

ORIENTATION 

PLV 

90. 

NUNBER 

9 

ORIENTATION 

PLV 

90. 

NUNBER 

10 

ORIENTATION 

PLV 
0 . 

NUNBER 

1 1 

ORIENTATION 

PLV 
0 . 

NUNBER 

12 

ORIENTATION 

PLV 

90. 

NUNBER 

13 

ORIENTATION 

PLV 

90. 

NUNBER 

14 

ORIENTATION 

PLV 
0 . 

NUNBER 

15 

ORIENTATION 

PLV 
0 . 

NUNBER 

16 

ORIENTATION 


SHOULD EACH PLV HAUE THE SAME SET OF ENGINEERING CONSTANTS? (V/N) 
V 


LONGITUDINAL VOUNGS NODULUS EX - 
19.5E6 

TRANSUEASE VOUNGS NOOULUS EV - 
1 . 32E6 

LONGITUDINAL TO TRANSUERSE POISSON RATIO NUXV - 
. 30 

SHEAR NOOULUS GXV - 
I . 0 I E6 


SHEAR NOOULUS GXZ 
t . 0 1 E6 
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SHEAR MODULUS GVZ - 
0 . 50E6 

LONG. THERMAL (HVGRO) COEFF. OF EXP. RLPX(I) - 

0.50E-6 

TRRH. THERMAL (HVGRO) COEFF. OF EXP. RLPV(I) - 

18.E-6 


**********»***«***»***«**«**«*«*«»«**««*«***«««»**«***»«»**«*******»»*** 


THE FOLLOUING INPUT DESCRIBES THE LOAOING 
CONDITIONS ON THE PLATE: 

TEMPERATURE (HVGRO) DIFFERENCE FROM REF. DELTA T • 
-180. 


TRANSUERSE PRESSURE LOAOING DELTA P - 


3 . 


USE CONTACT LAU? (V/N) 
V 


CONTACT LAU COEFF CONI? 

I . E6 

CRITICAL STRAIN ENERGY RELEASE PER UNIT AREA - 

0.2 

RELATIUE GROUTH DIRECTION OA/DB - 
0 . 


THE RELATIUE LOAD MAGNITUDES IN THE PLATE 
COORDINATE SYSTEM (1-PAIME, 2-PRIME). 

(FOR EXAMPLE, BN1 - t, BN2 - 0, BN6 - 0 
IS A SINGLE LOAO APPLIED IN THE I -PR I ME 
DIRECTION. ) 

BN f - 

1 .0 

BN2 - 
0 . 

BN6 - 

0 . 
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AT HOU NANV LOCATIONS IN THE OELAN I NAT I ON 00 VOU 
UISH STAESS/STAAIN CALCULATIONS? 

NUST - 

1 

INPUT 1 PAIAS OF COOAD I NATES : 

L I STX( 1) - 
0 . 

L I STV( 1) - 
0 . 

UOULD VOU LIKE THE STRESSES/STRAINS TO BE CALCULATED 
AT THE TOP (T), NIOOLE (H), OR BOTTON (B) OF EACH PLV? 

THE THROUGH-THICKNESS COORDINATE ORIGINATES UITH THE FIRST 
PLV OF THE SUBLAN I NATE AND PROCEEDS UNTIL THE LAST PLV OF THE 
PLATE. THUS, THE BOTTON OF THE FIRST PLV IS THE SUBLAHINATE 
OUTER SURFACE, AND SO ON. 


THE FOLLOUING INPUT PRESCRIBES UHICH CODE 
ANALYSIS OPTIONS UILL BE RUN: 

NONLINEAR LOAD-STRAIN HISTORY? (V/N) 

V 

OUTPUT STRAIN FILE HAHE? 

SANPLE. STRAIN 

OUTPUT STRESS FILE NANE? 

SANPLE. STRESS 

NAXINUN LOAD FOR POSTBUCKLING PLOT? 

-3000. 

NUNBER OF LOAD INCRENENTS FOR PLOT? 

5 

CALCULATE POSTBUCKLING GROUTH LOAD? (Y/N) 

V 


CALCULATE LINEAR 
V 

BUCKLING 

LOAD? 

(V/N) 

SHOULD 

V 

THE INPUT 

DATA BE 

URITTEN 

TO A FILE FOR FUTURE USE? (V/N) 

OUTPUT 

FILE NANE? 




SANPLE. 1 

INPUT 
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************************************************************************ 
INITIAL 2 COORDINATES 
PLY 2 

1 0 . OOe+OO 

2 0 . 56e-02 

3 0. 1 la-01 

4 0 . 1 7e-01 

5 0.22e-01 

6 0 . 28e-01 

7 0 . 33e-01 

6 0 . 39«-01 

9 0 . 44e-01 

10 0 . 50e-01 

11 0 . 56e— 0 1 

12 0 . 61 e-01 

13 0 . 67«-01 

14 0 . 720-01 

15 0 . 78«-01 

16 0 . 83e-01 

17 0.89e-01 


PLANE STRESS REDUCED STIFFNESSES FOR EACH PLV 


PLY 

QXX 

QVV 

QXY 

OSS 

QX2X2 

0Y2V2 

1 

0 . 20e+08 

0. 13e*07 

0 . 40e + 06 

0. 10e+07 

0. 10e+07 

0 . 50«+06 

2 

0 . 20e+08 

0 . 1 3e+07 

0 . 40e+06 

0 . 1 Oe+07 

0. 10e*07 

0 . 50e + 06 

3 

0 . 20«+08 

0.13e*07 

0.40«+06 

0 . 1 0e*07 

0 . 1 0a+07 

0 . 50e+06 

4 

0 . 20e+08 

0 . 1 3e*07 

0 . 40«+06 

0 . 1 0«*07 

0 . 1 0e*07 

0 . 50«+06 

5 

0 . 20e+08 

0. 13e+07 

0.40®*06 

0 . 1 0e+07 

0. 1 0e*07 

0 . 50e*06 

6 

0 . 20«+08 

0 . 1 3«+07 

0.40e>06 

0 . 1 0e+07 

0. 10e*07 

0 . 50e+06 

7 

0 . 20e+08 

0 . 1 3e*07 

0 . 40e+06 

0. 10®+07 

0. 10®*07 

0 . 50e*06 

8 

0 . 20e+08 

0 . 1 3«+0? 

0.40«+06 

0 . 1 0«+07 

0 . 1 0e+07 

0 . 50e>06 

9 

0.20e*Q8 

0. 13«*07 

0.40e*06 

0. 10«*07 

0. 10«+07 

0 . 50e*06 

10 

0.20**08 

0. 13«*07 

0 . 40«+06 

0. 10e*0? 

0. 10e*07 

0 . 50e+06 

1 1 

0 . 20e+08 

0 . 1 3®+07 

0 . 40«+06 

0 . 1 0a+07 

0 . 1 0«+07 

0 . 50«+06 

12 

0.20«+08 

0.13®*07 

0.40«*06 

0. 10«*07 

0 . 1 0®*07 

0.50**06 

13 

0.20«+08 

0 . 1 3«+0? 

0. 40®*06 

0. 10e*07 

0 . 1 0«+07 

0 . 50«+06 

14 

0.20#*08 

0. 13®*0? 

0 . 40e + 06 

0. 10e+07 

0 . 1 0e+07 

0.50e+06 

15 

0.20«+08 

0. 13«+07 

0.40e+06 

0. 10e*07 

0 . 1 0e+07 

0 . 50e+06 

16 

0 . 20«+08 

0.13«*0? 

0 . 40e+06 

0. 10e*07 

0 . 1 0e+07 

0 . 50e+06 


OFF-AX I S REDUCED STIFFNESSES IN THE PLATE COORDINATES 
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PLV 

on 

022 

012 

066 

Q16 

Q26 

1 

0 . 20e+08 

0. 13e+07 

0 . 40e+06 

0. 10e+07 

0. OOe+OO 

0. OOe + OO 

2 

0 . 20e+08 

0. 13e+07 

0 . 40e+06 

0 . 1 0e+07 

0. OOe+OO 

0. OOe+OO 

3 

0. 13e+07 

0 . 20e+08 

0 . 40e+06 

0. 1 0e+07 

0. OOe+OO 

0. OOe+OO 

4 

0. 13e+07 

0 . 20e+08 

0 . 40e+06 

0. 10e+07 

0. OOe+OO 

0. OOe+OO 

5 

0 . 20e+08 

0. 13e+07 

0 . 40e+06 

0 . 1 0e+07 

0. OOe+OO 

0. OOe + OO 

6 

0 . 20e+08 

0 . 1 3e+07 

0 . 40e+06 

0 . 1 0e+07 

0. OOe+OO 

0 . OOe+OO 

7 

0. 13e+07 

0 . 20e+08 

0 . 40e+06 

0. IOe+07 

0. OOe+OO 

0. OOe+OO 

8 

0. 13e+07 

0 . 20e+08 

0 . 40e+06 

0. 10e+07 

0. OOe+OO 

0. OOe+OO 

9 

0 . 1 3e+07 

0 . 20e+08 

0 . 40e+06 

0 . 1 0e+07 

0. OOe+OO 

0. OOe+OO 

10 

0 . 1 3e+07 

0 . 20e+Q8 

0 . 40e+06 

0. 10e+07 

0. OOe+OO 

0. OOe+OO 

1 1 

0 . 20e+08 

0. 13e+07 

0 . 40e+06 

0. 10e+07 

0. OOe+OO 

0. OOe+OO 

12 

0 . 20e+08 

0 . 1 3e+07 

0.40e+06 

0 . 1 0e+07 

0. OOe+OO 

0. OOe+OO 

13 

0. 13e+07 

0 . 20e+08 

0 . 40e+06 

0 . 1 0e+07 

0. OOe+OO 

0. OOe+OO 

14 

0 . 1 3e+07 

0 . 20e+08 

0 . 40e+06 

0 . IOe+07 

0. OOe+OO 

0. OOe+OO 

15 

0 . 20e+08 

0 . 1 3e+07 

0 . 40e+06 

0 . 1 0e+07 

0. OOe+OO 

0. OOe+OO 

16 

0 . 20e+08 

0 . 1 3e+07 

0 . 40e+06 

0 . 1 0e+07 

0. OOe+OO 

0. OOe+OO 


PLV 

044 

Q45 

Q55 

1 

0 . 50e+06 

0. OOe+OO 

0. 10e+07 

2 

0.50e+06 

0. OOe+OO 

0. 10e+07 

3 

0. 10e+07 

0 . OOe+OO 

0 . 50e+06 

4 

0. 10e+07 

0. OOe+OO 

0 . 50e+06 

5 

0 . 50e+06 

0. OOe+OO 

0 . 1 0e+07 

6 

0 . 50e+06 

0. OOe+OO 

0 . 1 0e+07 

7 

0 . 1 0e+07 

0. OOe+OO 

0 . 50e+06 

8 

0. 10e+07 

0. OOe+OO 

0.50e+06 

9 

0 . 1 0e+07 

0. OOe+OO 

0 . 50e+06 

10 

0. 10e+07 

0. OOe+OO 

0 . 50e+06 

11 

0 . 50e+06 

0. OOe+OO 

0 . 1 0e+07 

12 

0 . 50e+06 

0. OOe+OO 

0 . 1 0e+07 

13 

0. 10e+07 

0. OOe+OO 

0.50e+06 

14 

0. 10e+07 

0. OOe+OO 

0.50e+Q6 

15 

0 . 50e+06 

0. OOe+OO 

0. 10e+07 

16 

0.50e+06 

0. OOe+OO 

0. 10e+07 


OFF-AXIS REOUCEO STIFFNESSES IN THE SUBLANINATE COORDINATES 


PLV Qt 1 Q22 Q 1 2 066 


Q16 


026 


1 

2 

3 

4 


0 . 20e+08 
0. 20e+08 
0. 13e+07 
0 . 1 3e+07 


0 . 1 3e+07 
0.13e+07 
0 . 20e+08 
0 . 20e+08 


0 . 40e+06 
0 . 40e+06 
0 . 40e+06 
0.40e+06 


0. 10e+0? 
0. 10e+07 
0. 10e+07 
0. 10e+0? 


0 . OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


PLV 044 045 055 


1 

2 

3 

4 


0.50e+06 
0 . 50e+06 
0. IOe+07 
0. 10e+07 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0 . 1 0e+07 
0. IOe+07 
0 . 50e+06 
0 . 50e+06 



fin 1 NOTE : 

STIFFNESSES 

FOR THE 

PLATE 



mi 

022 

m 2 

A66 

A16 

A26 

0 . 93®+06 

0 . 93*406 

0.35e405 

0.90*405 

0.00*400 

0.000*00 

A44 

A45 

R55 




0 . 67*405 

0 . 00*400 

0.67*405 




B11 

022 

012 

B66 

016 

026 

0 . 73*-1 1 

0 . 00e+00 

-. 680-12 

0 . 45#-1 2 

0.00*400 

0.00*400 

044 

045 

055 




0 . 39*402 

0.00*4-00 

0.50*402 




F1 1 

F22 

FI 2 

F66 

F16 

F26 

0.11 *+01 

0 . 36*4-00 

0 . 28 e— 0 1 

0.70*-01 

0.00*400 

0.00*400 

F44 

F45 

F55 




0 . 42a-01 

O.OOe+OO 

0 . 63e-01 




H1 1 

H22 . 

H12 

H66 

HI 6 

H26 

0 . 1 7*-02 

0 . 36«-03 

0 . 39e-04 

0.990-04 

0.00*400 

0.00*400 

sn i hate ; 

STIFFNESSES 

FOR THE 

SUBLRIIINRTE 


m i 

022 

AI2 

A66 

fll 6 

826 

0 . 23e+06 

0.23*406 

0.89*404 

0.22*405 

0.00*400 

0.00*400 

(144 

045 

R55 




0 . 1 7*+05 

0.00*400 

0.17*405 




B11 

022 

012 

066 

016 

B26 

-.11 *+04 

0.11*404 

—.11 •- 1 2 

- . 2 8* — 1 2 

0 . 00*400 

0.00*400 


044 045 055 


0 . 69e+00 0 . 00e+00 0.69*400 
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0 . 71 e-03 0 . 71 e-03 0.27e-04 0.69e-04 0,00e*00 0.00e+00 

F44 F45 F55 

0.51 e-04 0 . 00e+00 0.51e-04 

Hit H22 HI 2 H66 H16 H26 

0 . 63e-07 0 . 63e-07 0.24e-08 0.61e-08 0.00e*00 0.00e+00 

**********************»*****t******************************************* 

LINEAR BUCKLING LORO CALCULATION: 

THE LINEAR BUCKLING LOAD IS -0. 67792e*03 

********************* *************************************************** 

GROUTH LOAD CALCULATION: 


UARNING: 

CONTACT 

BETUEEH 

SUBLAtllNATE 

AND 

PLATE 

UARNING: 

CONTACT 

BETUEEN 

SUBLAfl 1 NATE 

ANO 

PLATE 

UARNING: 

CONTACT 

BETUEEH 

SUBLAtllNATE 

AND 

PLATE 

UARNING: 

CONTACT 

BETUEEN 

SUBLAtllNATE 

AND 

PLATE 

UARNING: 

CONTACT 

BETUEEN 

SUBLAN 1 NATE 

ANO 

PLATE 

UARNING: 

CONTACT 

BETUEEN 

SUBLAtllNATE 

ANO 

PLATE 

UARNING: 

CONTACT 

BETUEEN 

SUBLAN 1 NATE 

AND 

PLATE 


THE GROUTH LOAD IS -0.17666e+04 ITER 


5 




STRESSES AND STRAINS AT (X,V) - ( 
FOR LOAD N - -0.17666e+04 


0 . OOOe+OO, 


0 . OOOe+OO ) 
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OFF-AXIS STRESSES IN THE SUBLAN I NATE FRANE 


PLV SI 


S2 S3 S4 S5 


1 

2 

3 

4 


17e+05 
- . 38e+05 
0 . 46e+03 
13e+04 


0 . 73e+04 
0.57e+04 
0. 18e+05 
0.22e+03 


0 . OOe+OO 
0 . OOe+OO 
0 . OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0 . OOe+OO 
0. OOe+OO 


OFF-AXIS STRAINS IN THE SUBLAN I NATE FRANE 


PLV El 


E2 


E3 


E4 


ES 


1 - . 1 0e-02 0.58e-02 0. OOe+OO 0. OOe+OO 0. OOe+OO 

2 - . 20e-02 0 . 49e-02 0. OOe+OO 0. OOe+OO 0. OOe+OO 

3 0 . 76e-04 0.90e-03 0. OOe+OO 0. OOe+OO 0. OOe+OO 

4 - . 96e-03 0.3U-04 0. OOe+OO 0. OOe+OO 0. OOe+OO 


ON-AXIS STRESSES IN THE PLV FRANE 

PLV SI S2 S3 S4 S5 


1 

2 

3 

4 


- . 1 7e+05 
- . 38e+05 
0 , 1 8e+05 
0 . 22e+03 


0 . 73e+04 
0 . 57e+04 
0 . 46e+03 
-. !3e+04 


0. OOe+OO 
0. OOe+OO 
0 . OOe+OO 
0 . OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


ON-AXIS STRAINS IN THE PLV FRAME 

PLV El E2 E3 E4 E5 


1 

2 

3 

4 


- . 1 0e-02 
- . 20e-02 
0 . 90e-03 
0.3le-04 


0 . S8e-02 
0 . 49e-02 
0 . 76e-04 
- . 96e-03 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0 . OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


OFF-AXIS STRESSES IN THE PLATE FRANE 

PLV St S2 S3 S4 S5 


1 

2 

3 

4 


17e+05 
- . 38e+05 
0 . 46e+03 
-. !3e+04 


0 . 73e+04 
0 . 57e+04 
0 . 1 8e+05 
0 . 22e+03 


0. OOe+OO 
0. OOe+OO 
0 .OOe+OO 
0. OOe+OO 


0 . OOe+OO 
0 . OOe+OO 
0. OOe+OO 
0 . OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


S6 

0 . OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


E6 

0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


S6 

0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


E6 

0. OOe+OO 
0. OOe+OO 
0 . OOe+OO 
0. OOe+OO 


S6 

0 . OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


OFF-AX I S MECHANICAL STRAINS IN THE PLATE FRAME 
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PLV El 


E2 E3 


E4 E5 


E6 


1 

2 

3 

4 


- . 75e-03 

0 . 29e-02 

0. OOe+OO 

0. OOe+OO 

0 . OOe+OO 

0 . OOe+OO 

-. 18e-02 

0 . 20e-02 

0. OOe+OO 

0. OOe+OO 

0. OOe+OO 

0. 

OOe+OO 

- . 28e-02 

0. 1 le-02 

0. OOe+OO 

0 . OOe+OO 

0 . OOe+OO 

0. 

OOe+OO 

- . 39e-02 

0 . 28e-03 

0. OOe+OO 

0. OOe+OO 

0. OOe+OO 

0. 

OOe+OO 


*******»»*••»****•*«*•*»***»»*»»*««»*•******««•****•*•******«****•****** 
NONLINERR LORD-STRAIN HISTORV CALCULATION: 


STRESSES AND STRAINS AT (X,V) - ( O.OOOe+OO, O.OOOe+OO) 
FOR LOAD N - O.OOOOOe+OO 

OFF-AXIS STRESSES IN THE SUBLANINATE FRANE 


PLV SI 


S2 S3 S4 SS 


S6 


1 

2 

3 

4 


- . 37e+04 
- ■ 38e+04 
0.38e+04 
0 . 38e+04 


0 . 38e+04 
0 . 38e+04 
- . 37e+04 
- . 38e+04 


0 . OOe+OO 
0. OOe+OO 
0 . OOe+OO 
0 . OOe+OO 


0 .OOe+OO 
0. OOe+OO 
0. OOe+OO 
0 . OOe+OO 


0. OOe+OO 
0. OOe+OO 
0 . OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


OFF-AXIS STRAINS IN THE SUBLANINATE FRANE 


PLV El E2 


E3 E4 


E5 


E6 


1 

2 

3 

4 


- . 25e-03 
- . 25e-03 
0 . 29e-02 
0 . 29e-02 


0 . 29e-02 
0 . 29e-02 
- . 25e-03 
-.25e-03 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0 . OOe+OO 


0. OOe+OO 
0. OOe+OO 
0 . OOe+OO 
0. OOe+OO 


ON-AXIS STRESSES IN THE PLV FRANE 


PLV SI 


S2 


S3 S4 S5 


S6 


1 

2 

3 

4 


-,37e+04 
- . 38e+04 
- . 37e+04 
-.38e+04 


0.38e+04 
0 . 38e+04 
0 . 38e+04 
0 . 38e+04 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0 . OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0 . OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


ON-AXIS STRAINS IN THE PLV FRANE 

PLV El E2 E3 E4 E5 E6 


1 

2 

3 

4 


- . 25e-03 
- . 25e-03 
- . 25e-03 
- . 25e-03 


0 . 29e-02 
0 . 29e-02 
0 . 29e-02 
0 . 29e-02 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0 . OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


OFF-AXIS STRESSES IN THE PLATE FRAHE 
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PLV SI 


S2 S3 


$4 


S5 


1 

2 

3 

4 


- . 37e+04 
-.38e+04 
0 . 38e+04 
0 . 38e+04 


0 . 38e+04 
0 .38e+04 
- . 37e+04 
- . 38e+04 


0 . OOe+OO 
0 .OOe+OO 
0 . OOe+OO 
0 .OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


OFF-RXIS flECHANICAL STRAINS IN THE PLATE FAAHE 


PLV El 


E2 


E3 E4 E5 


1 

2 

3 

4 


0. lle-06 
-.SOe-07 
- . 27e-06 
- .47e-06 


0 . 80e-06 
0 . 47e-06 
0 . I 3e-06 
- . 21 e-06 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


STRESSES AND STRAINS AT <X,Y) - ( 0.000e+00, 

FOR LOAD N - -0.60000e+03 


OFF-AXIS STRESSES IN THE SUBLANINATE FRANE 


PLV SI 


S2 S3 S4 SS 


1 

2 

3 

4 


!6e+05 
!6e+05 
0 . 29e+04 
0.29e+04 


0.3Se+04 
0.35e+04 
- . 35e+04 
-.35e+04 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


OFF-AXIS STRAINS IN THE SUBLANINATE FRANE 


PLV El E2 


E3 E4 ES 


1 

2 

3 

4 


- . 89e-03 
- . 90e-03 
0 . 23e-02 
0 . 23e-02 


0.29e-Q2 
0 . 29e-02 
-.23e-03 
-.23e-03 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


ON-AXIS STRESSES IN THE PLV FRANE 

PLV SI $2 S3 S4 SS 


1 

2 

3 

4 


-.!6e+05 0.35e+04 
- . I 6e+05 0.356+04 
- . 35e+04 0.29e+04 
- . 35e+04 0 . 29e+04 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 0. OOe+OO 
0. OOe+OO 0. OOe+OO 
0. OOe+OO 0. OOe+OO 
0. OOe+OO 0. OOe+OO 


S6 

0. OOe+OO 
0 . OOe+OO 
0. OOe+OO 
0. OOe+OO 


E6 

0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0.000e+00) 


S6 

0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


E6 

0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


S6 

0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 



ON-AX I S STRAINS IN THE PLV FRAME 
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PLV El 


E2 


E3 E4 E5 


1 

- . 89e-03 

0 

2 

- . 90e-03 

0 

3 

- . 23e-03 

0 

4 

- . 23e-03 

0 


29e-02 0 . OOe+OO 
29e-02 0. OOe+OO 
23e-02 0. OOe+OO 
23e-02 0. OOe+OO 


0. OOe+OO 0. OOe+OO 
0. OOe+OO 0. OOe+OO 
0. OOe+OO 0. OOe+OO 
0. OOe+OO 0. OOe+OO 


OFF-AXIS STRESSES IN THE PLATE FRAME 

PLV SI $2 S3 S4 $5 


1 

2 

3 

4 


16e+05 
- . 1 6e+05 
0 . 29e+04 
0 . 29e+04 


0 . 35e*04 
0 . 35e+04 
- . 35e+04 
- . 35e+04 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


OFF-AXIS MECHANICAL STRAINS IN THE PLATE FAANE 
PLV El E2 E3 E4 E5 


1 

2 

3 

4 


- . 64e-03 
- . 64e-03 
- . 65e-03 
- . 65e-03 


0 . 25e-04 
0 . 25e-04 
0 . 25e-04 
0 . 24e-04 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


STRESSES AND STRAINS AT <X,V) - ( 0.000e+00, 

FOR LOAD N - -0.12000e+04 


OFF-AXIS STRESSES IN THE SUBLAM I NATE FRANE 


PLV SI 


S2 S3 S4 SS 


1 

2 

3 

4 


- . 29e+05 
- . 29e+05 
0 . 21 e+04 
0 . 21 e+04 


0 . 33e+04 
0 . 33e+04 
- . 33e+04 
- . 33e+04 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


OFF-AXIS STRAINS IN THE SUBLAM I NATE FRAME 


PLV El E2 E3 


E4 


ES 


1 

2 

3 

4 


-. 15e-02 
-. 15e-02 
0. 16e-02 
0 . 1 6e-02 


0 . 29e-02 
0 . 29e-02 
- . 20e-03 
- . 20e-03 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


E6 

0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


S6 

0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


E6 

0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0 . 000e+00) 


S6 

0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


E6 

0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
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PLV SI S2 


S3 SI S5 


1 

2 

3 

4 


- . 29e+05 0 . 33e+04 
- . 29e+05 0 . 33e+04 
- . 33e+04 0 , 2 1 e+04 
-.33e+04 0 . 21 e+04 


0 .OOe+OO 
0 . OOe+OO 
0 . OOe+OO 
0 .OOe+OO 


0 . OOe+OO 
0 . OOe+OO 
0 . OOe+OO 
0 . OOe+OO 


0 , OOe+OO 
0. OOe+OO 
0 . OOe+OO 
0 . OOe+OO 


ON-AXIS STRAINS IN THE PLV FRANE 

PLV El E2 E3 E4 E5 


1 

2 

3 

4 


- . 1 5e-02 0 . 29e-02 
15e-02 0 . 29e-02 
- . 20e-03 0 . 1 6e-02 
- . 20e-03 0. 16e-02 


0 , OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 0. OOe+OO 
0. OOe+OO 0. OOe+OO 
0. OOe+OO 0. OOe+OO 
0. OOe+OO 0. OOe+OO 


OFF-AXIS STRESSES IN THE PLATE FRANE 

PLV SI S2 S3 S4 SS 


1 

2 

3 

4 


- . 29e+05 
- . 29e+05 
0.21e+04 
0.21e+04 


0 . 33e+04 
0 . 33e+04 
- . 33e+04 
-.33e+04 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


OFF-AXIS NECHANICAL STRAINS IN THE PLATE FRANE 
PLV El E2 E3 E4 ES 


1 

2 

3 

4 


- . 1 3e-02 
-.13e-02 
-. t3e-02 
-. !3e-02 


0.50e-04 
0 . 50e-04 
0.49e-Q4 
0.49e-04 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 0. OOe+OO 
0. OOe+OO 0. OOe+OO 
0. OOe+OO 0. OOe+OO 
0. OOe+OO 0. OOe+OO 


UARNING: CONTACT BETUEEN SUBLAfl I NATE AND PLATE 


STRESSES ANO STRAINS AT (X,V) - ( 0.000e+00, 

FOR LOAD N - -0.18000e+04 


OFF-AXIS STRESSES IN THE SUBLANIHATE FRANE 


PLV SI 


S2 S3 S4 


S5 


1 

2 

3 

4 


- . I 7e+05 
- . 30e+OS 
0.41 e+03 
-. !3e+04 


0.74e+04 
0. 58e+04 
0. IBe+05 
0.67e+03 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


S6 

0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


E6 

0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


S6 

0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


E6 

0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0.000e+00) 


S6 

0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 



OFF-AX I S STRAINS IN THE SUBLflfllNflTE FRAME 
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PLV El 


E2 E3 E4 ES 


1 

2 

3 

4 


- . 1 0e-02 
- . 21 e-02 
0 . 30e-04 
-. 10e-02 


0 . 58e-02 
0 . 50e-02 
0 . 93e-03 
0 . 55e-04 


0 .OOe+OO 
0 . OOe+OO 
0 . OOe+OO 
0. OOe+OO 


0. OOe+OO 
0 . OOe+OO 
0 . OOe+OO 
0. OOe+OO 


0. OOe+OO 
0 . OOe+OO 
0 . OOe+OO 
0. OOe+OO 


ON-RXIS STRESSES IN THE PLV FRAME 

PLV SI S2 S3 S4 S5 


1 

2 

3 

4 


17e+05 
- . 38e+05 
0. 18e+05 
0 . 67e+03 


0 . 74e+04 
0 . 58e+04 
0.41e+03 
- . 1 3e+04 


0. OOe+OO 
0 . OOe+OO 
0. OOe+OO 
0 . OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0 . OOe+OO 
0. OOe+OO 
0 . OOe+OO 
0. OOe+OO 


ON-AXIS STRAINS IN THE PLV FRAME 

PLV El E2 E3 E4 ES 


1 

-. 10e-02 

0 

2 

- . 21 e-02 

0 

3 

0 . 93e-03 

0 

4 

0 . 55e-04 

- 


58e-02 0. OOe+OO 
50e-02 0. OOe+OO 
30e-04 0. OOe+OO 
i0e-02 0. OOe+OO 


0. OOe+OO 0. OOe+OO 
0. OOe+OO 0. OOe+OO 
0. OOe+OO 0. OOe+OO 
0. OOe+OO 0. OOe+OO 


OFF-AXIS STRESSES IN THE PLATE FRAME 

PLV SI S2 S3 S4 S5 


1 

2 

3 

4 


-. 17e+05 
- . 38e+05 
0 . 41 e+03 
-. !3e+04 


0 . 74e+04 
0.58e+04 
0 . 1 8e+05 
0 . 67e+03 


0. OOe+OO 
0 .OOe+OO 
0 . OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


OFF-AXIS MECHANICAL STRAINS IN THE PLATE FRAME 


PLV El 


E2 E3 E4 


E5 


1 

2 

3 

4 


- . 75«-03 
- . !8e-02 
- . 29e-02 
- . 39e-02 


0 . 29e-02 
0 . 21 e-02 
0. f 2e-02 
0.3te-03 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 0. OOe+OO 
0. OOe+OO 0. OOe+OO 
0. OOe+OO 0. OOe+OO 
0. OOe+OO 0. OOe+OO 


UARNING: CONTACT BETUEEN SUBLAMINATE ANO PLATE 


E6 

0. OOe+OO 
0. OOe+OO 
0 . OOe+OO 
0 . OOe+OO 


S6 

0. OOe+OO 
O.OOe+QO 
0. OOe+OO 
0. OOe+OO 


E6 

0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


S6 

0 . OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


E6 

0. OOe+OO 
0 . OOe+OO 
0. OOe+OO 
0 . OOe+OO 


STRESSES AND STRAINS RT (X,V) - ( 0.000e+00, 0.000e+00) 

FOR LORO N - -0.24000e+04 

OFF-RXIS STRESSES IN THE SUBLRNINRTE FRRNE 


PLV SI 


S2 S3 $4 S5 


S6 


1 

2 

3 

4 


-. 19e+05 
- . 47e+05 
- . 38e+03 
- . 27e+04 


0 . 85e+04 
0 . 66e+04 
0 . 30e+05 
0. 92e+04 


0 . OOe+OO 
0. OOe+OO 
0 . OOe+OO 
0 . OOe+OO 


0. OOe+OO 
0 . OOe+OO 
0 . OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
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OFF-RXIS STRAINS IN THE SUBLRNINRTE FRRNE 


PLV El E2 


E3 E4 E5 


1 

2 

3 

4 


-.11 e-02 
- . 25e-02 
- . 74e-03 
- . 22e-02 


0.67e-02 
0 . 57e-02 
0. 15e-02 
0.51 e-03 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


ON-AXIS STRESSES IN THE PLV FRRNE 

PLV SI S2 S3 S4 $5 


1 

2 

3 

4 


-. 19e+05 
- . 47e+05 
0.30e+05 
0 . 92e+04 


0.85e+04 
0 . 66e*04 
- . 3Be+03 
-.27e+04 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


ON-AXIS STRAINS IN THE PLV FRRNE 

PLV El E2 E3 E4 E5 


1 

2 

3 

4 


-.11 e-02 
-.25e-02 
0. t5e-02 
0.51 e-03 


0.67e-02 
0 . 57e-02 
-.74e-03 
- . 22e-02 


0 .OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


OFF-RXIS STRESSES IN THE PLATE FRRNE 

PLV SI S2 S3 S4 S5 


1 

2 

3 

4 


- . 1 9e+05 
-.47e+05 
-.38e+03 
- .27e+04 


0. 85e+04 
0 . 66e+04 
0 . 30e+05 
0.92e+04 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 0. OOe+OO 
0. OOe+OO 0. OOe+OO 
0. OOe+OO 0. OOe+OO 
0. OOe+OO 0. OOe+OO 


E6 

0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


S6 

0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


E6 

0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


S6 

0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


OFF-ftXIS MECHANICAL STRAINS IN THE PLATE FAROE 
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PLV El 


E2 E3 E4 


E5 


E6 


1 

2 

3 

4 


- . 84e-03 
- . 22e-02 
- . 36e-02 
-.5le-02 


0 . 38e-02 
0 . 28e-02 
0 . 1 8e-02 
0 . 76e-03 


0 . OOe+OO 
0 , OOe+OO 
0 . OOe+OO 
0 . OOe+OO 


0. OOe+OO O.OO 0 +OO 
0. OOe+OO 0. OOe+OO 
0. OOe+OO 0. OOe+OO 
0. OOe+OO 0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


UARNIHG: CONTACT BETUEEH SUBLAtllHATE AND PLATE 


STRESSES AND STRAINS AT (X,V) - ( 0.000e+00, 

FOR LOAD N » -0.30000e+04 


OFF-AXIS STRESSES IN THE SUBLANINATE FRAME 


PLV SI 


S2 S3 S4 S5 


1 

2 

3 

4 


-.21e+05 
-.55e+05 
-.lle+04 
- . 38e+04 


0 . 93e+04 
0.72e+04 
0.40e+05 
0. IBe+OS 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


OFF-AXIS STRAINS IN THE SUBLAMINATE FRAME 


PLV El E2 


E3 E4 ES 


1 

2 

3 

4 


- . 1 2e-02 0.74e-02 0. OOe+OO 
- . 29e-02 0.63e-02 0. OOe+OO 
-. t5e-02 0.21e-02 0. OOe+OO 
- . 32e-02 0.99e-03 0. OOe+OO 


0. OOe+OO 0. OOe+OO 
0. OOe+OO 0. OOe+OO 
0. OOe+OO 0. OOe+OO 
0. OOe+OO 0. OOe+OO 


ON-AXIS STRESSES IN THE PLV FRAME 

PLV SI S2 S3 S4 S5 


1 

2 

3 

4 


- . 21 e+OS 
- . 55e+05 
0.40e+05 
O.I8e+05 


0.93e+04 

0.72e+04 

-.lle+04 

-.30e+04 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


ON-AXIS STRAINS IN THE PLV FRAME 

PLV El E2 E3 E4 E5 


1 

2 

3 

4 


- . 12e-02 
- . 29e-02 
0.2te-02 
0 . 99e-03 


0.74e-02 
0 . 63e-02 
-. ISe-02 
- . 32e-02 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


0.000e+00) 


$6 

0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


E6 

0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


S6 

0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 


E6 

0. OOe+OO 
0. OOe+OO 
0. OOe+OO 
0. OOe+OO 












OFF 

-AXIS STRESSES IN THE PLATE 

FRANE 
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• 

PLV 

St 

S2 

S3 

S4 

SS 

S6 


1 

- . 21 e+05 

0.93e+04 

0 . OOe+OO 

0 . OOe+OO 

0. OOe+OO 

0 . OOe+OO 



2 

- . 55e+05 

0 . 72e+04 

0. OOe+OO 

0 . OOe+OO 

0. OOe+OO 

0 . OOe+OO 



3 

1 1e+04 

0 . 40e+05 

0 . OOe+OO 

0 . OOe+OO 

0. OOe+OO 

0 . OOe+OO 



4 

-.38e+04 

0. 18e+05 

0 . 00e+00 

0. OOe+OO 

0. OOe+OO 

0 . OOe+OO 



OFF 

-AXIS NECHAN 1 CAL STRAINS IN 

THE PLATE 

FRANE 




PLV 

El 

E2 

E3 

E4 

ES 

E6 



1 

- . 96e-03 

0 . 45e-02 

0 . OOe+OO 

0. OOe+OO 

0. OOe+OO 

0. OOe+OO 



2 

- . 27e-02 

0 . 34e-02 

0 . OOe+OO 

0. OOe+OO 

0. OOe+OO 

0. OOe+OO 



3 

- . 44e-02 

0 . 23e-02 

0 . OOe+OO 

0. OOe+OO 

0. OOe+OO 

0. OOe+OO 



4 

* * * * i 

- . 61 e-02 

0 . 1 2e-02 

0. OOe+OO 

0. OOe+OO 

0. OOe+OO 

0. OOe+OO 



THE LINEAR BUCKLING LOAD OF THE DELANINATION IS -0. 67792e+03 


THE GROUTH LOAD OF THE DELANINATION IS 


-0. 17666e+04 




Appendix H 


Total Potential Energy Change of an Isotropic Plate 


The change in the total potential energy G of an isotropic, circular plate sub- 
jected to a transverse pressure q 0 with respect to a change in the area A = ira 2 of 
the plate is 


G = 


<m i 


da 2xa 

where the total potential energy II is calculated from 


(H.l) 


2tt a 


2ir 4 a 

U = ^J J J (a r e r + (r t e t )r dr dz dS — J J q 0 wr dr d$ 


(H. 2 ) 


0 0 00 

The plate thickness is ft, the radius a, and (r, $, z) are cylindrical coordinates. The 
transverse displacement of the plate is w. The radial and circumferential stresses, 
respectively, are calculated as 


Nr z .cPw 1 dw. 

<Tr ~~h~ 12D h?( dr 2 +U r dr^ 


{H. 3 ) 


N t z . cPw 1 dw. 

at ~T~ 12P ft3 1 Ir 2 + r dr ^ 


(HA) 
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where N r and JV t axe the radial and circumferential stress resultants and the bending 


stiffness D of the plate is 


D = 


Eh 3 

12(1 - i/ 2 ) 


(H.i) 


E is Young’s modulus, and u is Poisson’s ratio. 


The radial and circumferential 


strains, respectively, are 


Cr 





(*•«) 


__ u r z dw 

r r dr 

where u r is the radial displacement of the plate. The circumferential stress resultant 
N t is related to the radial stress resultant by 

N, = ±(rK) ( B.S ) 

and the radial displacement u r is related to the radial stress resultant by 

r dN 

Ur = jh {r ^r + ~ u)Nr) (jy - 9) 

A perturbation solution is developed by Chia [41] by expanding the transverse 
displacement u>, the pressure q 0 , and the radial stress resultant N r in terms of the 
displacement at the center of the plate w 0 . The transverse displacement is expanded 
as 



i/ W 0 W 0 3 

w = h ^ Wl ~h + ^ 


(/f.10) 
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where 


wi = £ 2 


(H- 11) 


= 55o (1 - ,/!) ^ <1 ■ (>< 8 \-T + 23{ + 8(2 + 2{S) (H12) 


and 


e-i--? 


a* 


(H. 13) 


The radial stress resultant is expanded in terms of the center displacement as 


AT Eh? . . Wo X 2 . / Wo ^4 v 

K = —Mfr + *«< — )) 


(*•14) 


where 


S2 = k 


6 1 — nu 




(*• 15) 


and 


(1-^/160-104*/ 80 - 52 v , 2 * 

, 4 ““roeo"V l-^ + "TrTT^ + ^ + * } 

- 501 ^ 249 */ 4 _ x 5 _ 39 ^6 _ h t \ 
l-u ) 


(IT. 16) 


The transverse pressure load is expanded in terms of the center displacement 


as 


,w 


*• ~ 3<r^ ( T + So (I + v)( ™ - nv ^T )3) 


{H. 17) 


For a Poisson’s ratio v = 0.3, the change in the total potential energy (Eq. H.l) 
was evaluated as a function of the normalized center displacement as 


Ga 4 


= 2.930(^) 2 +l,586(^) 4 +0.127S(^) 5 +0.003048(^) , +8.0064£-5(^) 

(H. 18) 


, w 


o \10 
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Similarly, the transverse pressure (Eq. H.17) was evaluated as 


Eh 4 


= 5.860(^) + 3.657(^) 3 
h h 


(ff.19) 


These equations were then parametrically evaluated to determine the relationship 
between the applied transverse pressure and the change in total potential energy. 



Appendix I 


Ultrasonic Nondestructive Examination 


Every specimen in the experimental portion of this dissertation was ultrason- 
ically examined before and after compression testing. The ultrasonic scanning, 
commonly known as C-scanning, was performed in the Structures and Composites 
Laboratory on equipment built and programmed by the author. The C-scans pro- 
vided a planar map of each specimen showing the lateral extent of delamination and, 
in particular, the depth in number of plies of the delamination at every point. The 
data were used before testing to precisely locate the teflon implants to apply strain 
gauges to the specimen surface, and to map the extent of delamination growth after 
testing. 

The C-scan equipment consists of: (a) an ultrasonic flaw detector (Kraut- 
kramer Branson USL 48), (b) an immersion tank and specimen positioning fix- 
tures, (c) a bridge with stepper motors to drive the ultrasonic transducers back 
and forth over the specimen (Trienco Model 705), and (d) a computer to perform 
data acquisition and control functions as well as to display the color output (IBM 
PC/AT with IBM data acquisition card). The C-scan was operated by a FOR- 
TRAN computer program that controlled the movement of the transducer bridge, 
collected the data, and converted the data into a graphical display. 

The C-scan operation is based on generating a pulse of ultrasonic sound by a 
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transducer. The pulse travels through a coupling medium (water) to the specimen. 
At every interface between two dissimilar media, part of the signal will be reflected 
and part transmitted. Thus, there will be reflected signals from the top and bottom 
surfaces of the specimen as well as from any delaminated surfaces in between. In 
the pulse/echo method, the first signal returning from the top surface is used as 
a trigger and the time for subsequent signals to arrive is measured. Knowing the 
speed of sound in the material, the time of flight measurements axe converted to 
thicknesses. The thicknesses axe finally displayed as depths to the delamination at 
that point or, if there is no delamination, as the overall thickness of the specimen. 



Appendix J 


Uncertainty Analysis 


§J.l Experimental Uncertainty 


Uncertainty in the experimental data occurs due to variations in specimen fab- 
rication, preparation, strain gauging, testing, data acquisition, and data reduction. 
The purpose of this analysis is to estimate the uncertainty in a measured load 
associated with a given value of strain. Each experiment had four strain gauges 
mounted away from the delamination whose purpose was to measure the far field 
strain (gauges 2, 3, 4, and 9, Figure 6-2, Chapter 6). The response of these gauges 
should nominally be the same for a given experiment series, and thus may serve as 
replicate strain readings. For example, Figure J-l shows the measured load versus 
strain from gauge 3 for each of the four experiments in Series 5. 

The method used to calculate the experimental uncertainty is to first fit a linear 
least squares regression line to the data, and then to estimate the experimental 
uncertainty from the differences, or residuals, between the regression line and each 
data point [55, 56]. The estimate of the data experimental uncertainty, crj, is 
calculated from 


<?i 


2 


NM - naitii 

n — 1 


(j.i) 
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where a* is one standard deviation, is the applied load per unit specimen width 
corresponding to a single data point, Cj the associated strain, and n the number of 
data points (repeated subscripts imply summation). The estimated experimental 
uncertainties for gauges 2, 3, and 4 from each test series are summarized in Table J- 
1. Gauge 9 was not included because it was transversely oriented, and therefore 
substantively different from gauge 1, the gauge of primary interest in the delamina- 
tion studies. The estimated uncertainty for Test Series 4 gauge 3 is very large. One 
gauge from this series was clearly different from the others, indicating a systematic 
and not random variation. 



Figure J-l Load versus strain from each gauge 3 of Experiment Series 5. 
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Gauge 2 

Gauge 3 

Gauge 4 

Test Series 4 

43.7 

604. 

52.5 

Test Series 5 

57.6 

63.8 

68.5 

Test Series 6 

150. 

126. 

142. 


§J.2 Prediction Uncertainty 

Random uncertainty in the model predictions is due to errors in the input data 
propagating through the code. The response surface method may be used to es- 
timate this uncertainty [57, 58, 59]. The model predictions N are evaluated for 
different combinations of perturbations to the input data Xi, termed the experi- 
mental design, where the perturbations are plus or minus one standard deviation 
<t i of the input variable about its nominal mean value m. The model responses so 
generated are used to fit a truncated Taylor’s series expansion in the input vari- 
ables, which is then used to estimate the prediction uncertainty. The Taylor’s series 
expansion of the model is 

dN 

N(xi) = N(tii) + —(Xi -/*) + ••• (.7.2) 

The model prediction uncertainty as for a calculated load N is estimated to first 
order from 

o ’ jv 2 = (^~ <T «) 2 (* 7 - 3 ) 

where the Oi are one standard deviations of the input variables. 

The uncertainty in the model prediction was estimated for the particular case 
of Experiment 6-2. Table J-2 lists the means and standard deviations for fourteen 
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input variables considered in the analysis. The experimental design was a 2 14-10 
fractional factorial design, meaning that 2 14 cases would have to be run to include 
every possible combination of plus and minus factors, but that only a 2 10 fraction of 
the full design was run (16 cases). In this case only the linear terms of the Taylor’s 
series expansion could be estimated. 

The model prediction uncertainties were estimated for four different loads: (a) 
the linear buckling load, (b) the nonlinear buckling load, (c) the growth load, and 
(d) the load at 1000 microstrain (postbuckling regime). The mean values of the 
model predictions and the associated one standard deviations are listed in Table J- 
3. The transverse pressure was specifically not included in the uncertainty analysis 
even though the model is known to be sensitive to it since there was no way of 
estimating the uncertainty in it. For a graphical sense of the model prediction 
uncertainty, Figure J-2 shows the load versus strain responses corresponding to the 
sixteen different cases run in the uncertainty analysis. 

The uncertainty in the growth load is dominated by the uncertainty in the 
critical strain energy release rate which, as discussed in Chapter 7, is not well char- 
acterized for the material used in these experiments. By contrast, the estimated 
prediction uncertainty in the growth load without a contribution from the critical 
strain energy release rate is 54.9 lbf/in. Similarly, the uncertainty in the postbuck- 
ling load is dominated by the strain gauge thickness uncertainty. Without this 
contribution the uncertainty in the load at 1000 microstrain is 177. lbf/in. 
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Table J-2 Input Variable Uncertainties in Experiment 6-2 


Variable ^ a 

Units 


1 

Contact law, K 

1.E6 

0.5E6 

in 3 

2 

Sublaminate major semi-axis, a 

1.000 

0.033 

in 

3 

Sublaminate minor semi-axis, b 

0.750 

0.033 

in 

4 

Sublaminate angle with respect to loads, 9 

30 

2 

degrees 

5 

Ply thickness, t 

5.56E-3 

0.093E-3 

in 

6 

Longitudinal Young’s modulus, E x 

19.5E6 

0.65E6 

psi 

7 

Transverse Young’s modulus, E y 

1.32E6 

0.04E6 

psi 

8 

Poisson’s ratio, u xy 

0.30 

0.01 

- 

9 

In-plane shear modulus, G xy 

1.01E6 

0.03E6 

psi 

10 

Longitudinal thermal coeff. of expansion, a x 

0.50E-6 

0.017E-6 

in ^ 

11 

Transverse thermal coeff. of expansion, a x 

18.0E-6 

0.6E-6 

in 

in— °F 

12 

Temperature change, AT 

-180 

20 

°F 

13 

Critical strain energy release rate, G c 

0.3 

0.05 

in—lbf 

in 2 

14 

Gauge thickness, t g 

0.003 

0.001 

in 


Table J-3 Prediction Uncertainty Analysis 



UN 

a n 

Units 

Linear buckling load, N l h 

619. 

58.1 

lbf/in 

Nonlinear buckhng load, Nb 

1012. 

58.6 

lbf/in 

Growth Load, N g 

1763. 

162. 

lbf/in 

Load at 1000 microstrain, 7V 10 oo 

1541. 

245. 

lbf/in 





Appendix J: Uncertainty Analysis 


136 



Figure J-2 Uncertainty analysis of Experiment 6-2 prediction. Load versus strain 
for the sixteen different combinations of input variables. 




